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1 .  Scope 


Numerical  models  have  been  an  effective  tool  in  the  prediction  of  many 
geophysical  systems.  Processes  within  the  earth’s  two  primary  physical  systems,  the 
ocean  and  atmosphere,  directly  impact  Naval  operations  on  the  mesoscale  and  the 
local  PBL  scale.  Consequently,  accurate  and  efficient  prediction  on  these  scales  is  a 
necessity.  Scientists  now  consider  the  ocean  and  atmosphere  as  separate  but  fully 
coupled,  two-way  interactive  fluids.  Predicting  the  behavior  of  either  fluid  depends  upon 
the  spatial  and  temporal  forcing  applied  by  the  other.  Thus,  a  single  numerical 
prediction  system  combining  an  oceanic  and  an  atmospheric  model  provides  more 
realistic  representation  of  these  two  geophysical  systems.  Additionally  in  the  past 
decade,  increased  computer  power  and  technological  advancements  have  improved 
computational  efficiency  allowing  larger  models,  with  higher  resolution,  multi-nested 
grids  and  complicated  physics,  to  be  developed  and  run  for  real-time  forecasting 
purposes.  When  used  in  a  research  mode,  the  models  also  provide  valuable  insight 
toward  understanding  complex  mesoscale  interactions.  To  fully  utilize  modern 
computer  resources  and  to  meet  the  growing  need  for  high  resolution,  coupled 
oceanic/atmospheric  forecasts,  a  new  model  has  been  developed  by  the  Naval 
Research  Laboratory:  The  Coupled  Ocean/Atmosphere  Mesoscale  Prediction  System 
(COAMPS). 

At  present,  COAMPS  consists  of  two  FORTRAN  programs:  an  analysis  program 
that  blends  observations  with  first  guess  fields  to  provide  initial  conditions,  and  a 
forecast  program  containing  a  nonhydrostatic,  quasi-compressible  atmospheric  model 
with  complete  physics  schemes  for  predicting  meso  and  micro  scales  of  motion  (i.e. 
time  scales  ranging  from  days  to  minutes  and  spatial  scales  ranging  from  thousands  of 
kilometers  to  a  few  meters).  This  system  can  be  run  using  idealized  or  real  data  initial 
conditions  with  up  to  seven  horizontally  nested  domains.  Because  the  ocean  forecast 
model  of  COAMPS  currently  is  under  development,  it  will  be  implemented  at  a  later 
date.  For  this  reason,  this  report  documents  only  the  atmospheric  model.  Moreover, 
the  focus  here  is  to  provide  a  brief  introduction  to  COAMPS,  along  with  a  quick 
reference  guide  for  setting  up  and  performing  a  COAMPS  model  run.  It  is  not  intended 
as  a  complete  description  of  the  model  physics  or  dynamics.  Additional  reading 
material  on  COAMPS  and  related  references  are  given  in  Sections.  An  overview  of  the 
COAMPS  directory  structures  and  the  code  is  given  in  Section  3.  in  Section  4.  the 
procedure  for  executing  a  model  run  is  described  for  either  idealized  or  real  data 
assimilation  process  studies.  This  section  is  beneficial  in  assisting  the  new  user  and 
also  as  a  reference  guide  for  identifying,  defining  and  modifying  relevant  model 
parameters.  Appendices  A-E  provide  sample  files  and  tables  that  list  and  define 
COAMPS  input  parameters. 
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3.  System  Overview 

To  effectively  utilize  all  the  features  of  COAMPS,  the  user  must  become  familiar 
with  two  aspects  of  the  system:  the  COAMPS  directory  structure  and  coding  framework, 
and  the  process  for  submitting  a  COAMPS  job.  While  the  latter  subject  is  covered  in 
Section  4,  here  a  brief  overview  of  the  complete  modeling  system  is  provided.  The 
initial  COAMPS  directory  tree  consists  of  several  subdirectories  that  contain  the 
COAMPS  source  code,  databases,  templates  for  setting  up  a  COAMPS  job, 
documentation,  and  utility  programs.  Discussion  of  this  directory  structure  is  followed 
by  an  outline  of  the  COAMPS  coding  which  is  made  up  of  an  analysis  and  forecast 
driver  program.  The  use  of  flow  chart  diagrams  and  subroutine  descriptions  provides 
an  overview  of  the  two  programs.  Both  the  COAMPS  directory  structure  and  coding 
framework  are  of  critical  importance  in  understanding  the  modeling  system,  working 
with  the  code  and  properly  executing  a  job. 

The  following  conventions  are  adopted  throughout  the  document:  all  COAMPS 
variables  (including  inputs,  array  names  and  other  model  parameters)  are  in  boldface 
type.  Appendix  C,  D  and  E  lists  and  defines  each  of  these  model  input  parameters  and 
is  to  be  used  as  a  reference  while  reading  this  document.  Filenames  (including 
subroutines,  executables,  script  files,  and  namelists)  are  in  italics.  Some  filenames 
contain  a  *$’  indicating  variable  substitution.  Before  describing  the  COAMPS  directory 
structure  and  coding,  a  few  important  general  features  are  introduced  below, 

3.1.  COAMPS  General  Features 

The  COAMPS  code  is  more  effectively  outlined  by  first  summarizing  a  few  details 
about  the  model.  Introduced  in  the  following  paragraphs  are  the  driver  programs, 
namelists,  model  domain,  data  assimilation  procedure  and  COAMPS  databases. 

3.1.1.  Driver  Programs  and  Namelists 

Because  COAMPS  is  a  large  system,  comprising  nearly  83,000  lines  of 
FORTRAN  code,  information  can  best  be  conveyed  by  flow  chart  diagrams  --  the  first  of 
which  is  shown  in  Figure  1  for  the  COAMPS  analysis  and  forecast  driver  programs, 
coamps_analysis.Fa.n6  coamps_forecasf.F  respectively.  The  driver  program’s  function 
is  to  read  in  the  user  defined  model  grid  dimensions  that  provide  array  allocation  space 
and  pointer  locations  for  all  of  COAMPS  spatially  dependent  variables.  These  variables 
consist  of  the  number  of  grid  nests  (nnest),  the  number  of  horizontal  grid  points  for 
each  nest  (ma{nn),  na(nn),  where  nn  Is  the  nest  number),  and  the  number  of  vertical 
levels  (kka).  Their  values  are  specified  by  reading  in  the  namelist  gridni.  Then 
subroutine  mema.Fls  called  to  compute  the  total  length,  nwords,  of  a  single  array  (a) 
that  allocates  array  space  for  all  the  analysis  fields  on  all  meshes.  In  addition  to 
computing  the  length  of  array  a,  mema.F also  assigns  the  pointer  location  within  that 
array  for  each  field  on  each  nest.  These  pointer  locations  are  denoted  by  the  variable 
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coamps_analysis.  f 

analysis  driver  program 


Figure  1 :  Flow  chart  of  the 
CO  AMPS  driver  programs 
coamps _analy sis. f  and 
coamps_forecast.f.  The  model 
domain  specifications  are  read 
in  through  the  gridnl  namelist 
and  the  pointers  and  array  space 
are  setup  in  subroutines  mema.f 
and  memm.f  before  calling  the 
main  analysis  and  forecast 
subroutines  coama.fand 
coamm.f. 


user  input 

gridnl  namelist 
model  domain  parameters 


coamps_forecast  f 

forecast  driver  program 


user  input 

gridnl  namelist 
model  domain  parameters 


name  preceded  by  an  ‘i’.  For  example,  a  3D  variable  named  ‘var’  that  is  typically 
dimensioned  var(ma(nn),na(nn),kka)  becomes  a(ivar(nn)),  where  Ivar  points  to  the 
first  position  in  array  a  that  contains  the  value  of  variable  var  for  nest  number  nn. 

Finally,  the  grid  dimensions,  and  the  array  allocation  and  pointer  information  are  passed 
to  the  primary  analysis  routine,  coama.F,  where  the  initial  forecast  fields  are  prepared. 
Once  these  fields  are  output  by  the  analysis  program,  the  forecast  driver  program, 
coamps_forecast.F,  performs  similar  steps  to  set  up  array  allocation  space  and  pointer 
information  that  is  passed  to  the  primary  forecast  routine,  coamm.F,  where  the 
integration  of  the  model  equations  generates  the  prognosis. 

In  addition  to  the  gridnl  namelist,  COAMPS  also  reads  in  user  input  through  two 
other  namelists:  coamnland  mvoini,  which  are  covered  in  greater  detail  in  several  of  the 
following  sections.  The  use  of  namelists  permits  coding  flexibility  by  allowing  the  user 
to  manipulate  COAMPS  features  at  execution  time  without  requiring  recompilation  of 
the  code.  Appendix  B  shows  an  example  of  these  namelist,  and  Appendices  C  and  D 
give  a  complete  list  of  all  the  possible  user-specified  inputs. 

3.1 .2.  Model  Domain  Structure 

The  COAMPS  system  can  be  run  using  up  to  seven  horizontally  nested  grids 
with  the  horizontal  resolution  for  each  inner  nest  increasing  by  three  times  that  of  the 
next  larger  nest.  Consequently,  as  shown  in  Figure  2,  every  third  grid  point  on  an  inner 
nest  is  coincident  with  a  grid  point  on  the  next  larger  nest,  referred  to  as  it’s  ‘parent 
mesh’.  In  the  triply  nested  example  shown  in  Fig.  2,  the  outer  nest,  denoted  the 
‘coarse’  mesh  (1),  has  horizontal  dimensions  determined  by  namelist  input  parameters 
ma(1),  na(1).  Correspondingly,  the  first  inner  nest  (2)  has  horizontal  dimensions  given 
by  ma(2),  na(2),  and  the  third  inner  nest  (3)  by  ma(3),  na(3}. 

On  any  given  nest,  COAMPS  uses  a  vertically  and  horizontally  staggered  grid, 
shown  schematically  in  Figures  3a  and  3b  respectively.  In  these  figures,  the 
momentum  components  (u,v,w)  are  shifted  one-half  grid  interval  from  the  location  of  the 
other  prognostic  variables  (e,9,7c,q’s).  The  vertical  staggering  requires  that  the  vertical 
velocities  (w)  be  defined  on  height  levels  computed  from  the  namelist  input  array 
dsigma  defined  in  Appendix  D.  This  array  contains  layer  thickness’  assigned  by  the 
user  in  the  coamn/ namelist.  In  Figure  3  and  in  the  COAMPS  code,  heights  computed 
from  these  layer  thickness’  are  represented  by  array  sigmwa.  All  other  prognostic 
variables  (u,v,e,0,7C,q’s)  are  located  halfway  between  two  sigmwa  levels  at  heights 
represented  by  array  sigmma.  Hereafter,  these  heights  are  denoted  more  generally  as 
‘sigma  levels’.  It  is  important  to  note  that  In  setting  up  COAMPS  arrays  the  vertical 
index  k  increases  from  the  top  down. 
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Coarse  mesh  (1) 


Figure  2;  Schematic  representation  of  three  horizontally  nested  domains.  The  outer  coarse  mesh  (1) 
contains  horizontal  grid  points  ma(l)=10  and  na(l)=8.  The  first  inner  nest  (2)  contains  horizontal 
grid  points  ma(2)=10  and  na(2)=10.  The  second  inner  nest  (3)  contains  horizontal  grid  points 
ma(3)=7  and  na(3)=7.  Note  that  every  third  inner  nest  grid  point  is  collocated  Avith  a  grid  point  on 
it’s  parent  nest,  denoted  by  the  open  circles.  Thus,  the  number  of  inner  nest  grid  points  - 1  must  be 
divisible  by  3  (e.g.  (ma(2)-l)/3=3). 


sigrawa(l) 
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Figure  3a:  Schematic  representation  of  the  vertical  grid  staggering.  Vertical  velocities  are  defined  on  sigmwa 
levels  computed  from  user-specified  height  increments  dsigma.  All  other  prognostic  variables  are  defined  on 
sigmma  levels  located  halfway  between  the  sigmwa  heights.  In  this  example,  (m,n,kk)=(5,5,8)  and  the 
vertical  plane  shown  is  fixed  at  j=3.  The  array  indices  for  one  row  and  column  of  variables  w(ij,k),  u(i,j,k), 

7l(i,j,k)  are  shown  for  reference,  where  11,  the  dimensionless  pressure,  is  representative  of  the  other  mass  variables. 


;r(1.3,8)  ii(l,3.8)  ;r(2,3,8)  «(2,3,8)  a-(3,3,8)  tt(3,3,8)  ;r(4.3.8)  u(4,3.8)  ^-(5,3,8)  Mf5.3.8) 


delx 


Figure  3b;  Schematic  representation  of  the  horizontal  grid  staggering. 
The  u-  and  v-momentum  components  are  shifted  one-half  grid  interval 
from  the  other  prognostic  variables.  The  delx  and  dely  are  the  x  and  y 
grid  intervals  specified  by  the  user  for  the  coarse  mesh.  In  this  example, 
(m,n,kk)=(5,5,8)  and  the  horizontal  plane  shown  is  fixed  at  k=8.  The 
array  indices  for  one  row  and  column  of  variables  u(i,j,k),  v(i,j,k),7l(i,j,k) 
areshown  for  reference,  where  7C,  the  dimensionless  pressure,  is 
representative  of  the  other  mass  variables. 
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3.1 .3.  Initial  Conditions  and  COAMPS  Databases 

To  begin  any  simulation,  the  initial  conditions  must  be  specified.  The  initial 
conditions,  also  referred  to  as  the  initial  forecast  fields,  are  comprised  of  several 
COAMPS  arrays  that  contain  values  for  the  surface  data,  basic  state  variables  and 
prognostic  fields  on  each  mesh.  The  values  are  given  by  the  analysis  program  from 
idealized  user  input  or  alternatively  from  real  data  obtained  from  surface  databases, 
standard  atmospheric  values,  first  guess  fields  and  MVOI  increments.  The  COAMPS 
arrays  that  correspond  to  these  initial  fields  are  given  in  Tables  1a-c. 


Table  la:  COAMPS  2D  surface  arrays.  Idealized  values  are  used 
when  the  appropriate  nsunelist  flag  is  set  to  0 ;  otherwise  the 
database  values  are  used. 


ARRAY 

DEFINITION 

NAMELIST 

NAMELIST 

CO AMPS 

FLAG 

VAR.  VALUE 
or  DATABASE 

UNITS 

xiand 

land/sea  table 

ilndflgsO-^ 

ilndflg=1-^ 

aingpct 

dsclim 

0=water; 

1  =land 

dimensionless 

zsfc 

terrain  height 

itopoflg=0-> 

itopoflg=1-4 

0.0 

dsclim, dsdted 

meters 

psfc 

surface  pressure 

ItopoflgsO-^ 

psfcO 

ImOSHH 

albedo 

albedo 

ialbflg=0-> 

ialbflg=1-> 

albdo 

dsclim, dsgiss 

0=absorb; 

1  =reflect 
dimensionless 

tsea 

surface  temperature 

iseaflg=0-> 

iseaflg=1-» 

seatmp 

Kelvin 

tsoil 

deep  soil  temperature 

— 

dsclim 

Kelvin 

gwet 

ground  wetness 

iwetflg=0-> 

iwetflg=1-^ 

sfcwet 

dsclim 

0=dry;  1=sat 
kg/kg 

zO 

surface  roughness. 

r'ffftW 

zrough 

dsclim,dsgiss 

meters 

Table  lb:  COAMPS  3D  basic  state  arrays.  Idealized  values  are 
used  when  namelist  variable  icase^iO;  Standard  atmospheric  values 
are  used  when  icase=0 . 


ARRAY 

DEFINITION 

NAMELIST 

FLAG 

NAMELIST  VAR 
or  VALUE 

ugeoa 

geostrophic  u-wind 
component 

icase?i0-^ 

icase=0-» 

usnd  &  umean 
ugeoa 

m/s  (iuvsnd=1) 
dir/spd  (luvsnd=2) 

vgeoa 

tbw/tbm 

geostrophic  v-wind 

component 

temperature 

icase^iO-^ 

icasesO-^ 

icase?H)-> 

vsnd  &  vmean 

vgeoa 

tsnd 

m/s  (iuvsnd=1) 
dir/spd  (iuvsnd=:2) 
Kelvin  (itsnd=1) 
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tbw/tbm 

qvbw  & 
qvbm 

temperature 

specific  humidity 

icase=0-^ 

lcase9i0-> 

icase=0-> 

standard  atmos 

qsnd 

computed  from 
standard  atmos 

Celsius  (itsnd=2) 
theta  (itsnd=3) 
dewpt  K  (iqsndsl) 
dewpt  C  (iqsnd=2) 
g/kg  (iqsnd=3) 
rel  hum  (iqsnd=4) 

prbw  & 
prbm 

pressure 

Icase^tO-^ 

icase=0-^ 

ipsnd 

standard  atmos 

Pascal:  kg/{ms^) 

exbw  & 
exbm 

dimensionless 

pressure 

computed 

Internally 

Pascal;  kg/(ms^) 

thbw  & 
thbm 

potential 

temperature 

computed 

internally 

Kelvin 

thvbw  & 
thvbm 

virtual  potential 
temperature 

computed 

internally 

Kelvin 

rbw  &  rbm 

density 

kg/m“ 

Table  Ic:  COAMPS  3D  prognostic  arrays.  Idealized  values  are 
used  when  namelist  variable  icase^^O;  First  guess  fields  (and  MVOI 
increments  when  iupd=2 )  are  used  when  lcase=0 . 


ARRAY 

DEFINITION 

Ul 

u  momentum 

v1 

V  momentum 

Wl 

vertical  velocity 

thi 

potential  temperature  (0) 

el 

turbulent  kinetic  energy  (TKE) 

qvl 

mixing  ratio:  water  vapor 

qc1/qi1/ 

qr1/qs1 

mixing  ratios: 
clouds/ice/rain/snow 

dimensionless  pressure 
perturbation  (tt) 

NAMELIST 


FLAG 


icase?iO-> 

icase=0->^ 


icasejiO-^ 

lcases=0-> 


icase9i0-» 

icase=0-» 


NAMELIST  VAR 


or  VALUE 


usnd  &  umean 

first  guess  field  & 
MVOI  increments 


vsnd  &  vmean 

first  guess  field  & 
MVOI  increments 


0.0 

first  guess  field 


thbm 

first  guess  field  & 
MVOI  increments 


1.0e-10 

first  guess  field 


qvbm 

first  guess  field  & 
MVOI  increments 


0.0 

first  guess  field 


0.0 

first  guess  field  & 
MVOI  increments 


UNITS 


Because  COAMPS  uses  the  leap-frog  integration  scheme,  three  time  levels  for 
each  of  the  3D  prognostic  fields  must  be  retained  during  each  iteration.  Thus,  the  array 
names  are  further  identified  by  a  time  level  number  1 , 2,  or  3: 

Table  2 :  COAMPS  array  naming  convention 


NUMBER 

TIME  LEVEL 

1 

previous  (t-At) 

2 

present  (t) 

3 

predicted  (t+At) 

Note  that  the  3D  arrays  shown  in  Table  1c  are  for  time  level  1  representing  values  from 
the  previous  iteration.  For  purposes  of  discussion,  the  arrays  are  denoted  here  and 
throughout  the  remainder  of  the  document  by  the  conventional  notation 
vartl(ma(1),na(1),kka).  In  the  code,  the  3D  prognostic  arrays  are  denoted  by 
a(ivarti(nn)),  where  ‘a’  is  the  forecast  program  array,  ‘i’  is  the  pointer  location  for 
variable  ‘var’  at  time  level  1:1’  on  nest  number  ‘nn’.  As  mentioned,  the  array  dimension 
parameters  (ma(nn),  na(nn),kka)  correspond  to  the  number  of  grid  points  for  each 
mesh  that  are  specified  by  the  user  in  the  gridnl  namelist.  In  general,  all  of  the  above 
arrays,  along  with  other  fields  and  variables  In  COAMPS,  are  in  MKS  units.  With  these 
preliminary  details  in  mind,  the  COAMPS  directory  structure  is  described  below. 

3.1.4.  First  Guess  Fields  and  Data  Assimilation 

Real  data  initial  conditions  require  a  set  of  first  guess  fields  that  represent  the 
initial  state  of  the  atmosphere.  These  fields  are  given  in  one  of  two  ways:  from  a 
NOGAPS  analysis  or  forecast,  defined  as  a  ‘NOGAPS  cold  start’,  or  from  a  previous 
COAMPS  forecast,  defined  as  a  ‘data  assimilation  update  cycle’. 

The  NOGAPS  cold  start  provides  global  fields  on  pressure  levels  and  at  1  ° 
horizontal  resolution  that  are  interpolated  to  the  model’s  horizontal  grid  points.  These 
fields  are  either  initial  conditions  to  the  forecast  program,  or  optionally,  a  first  guess  to 
the  analysis  program  where  a  multi-variate  optimum  interpolation  (MVOI)  scheme 
blends  in  observational  data. 

The  data  assimilation  update  cycle,  shown  schematically  in  Figure  4,  begins  with 
a  NOGAPS  cold  start  and  an  MVOI  analysis  to  produce  COAMPS  forecast  fields  valid 
one  or  two  days  prior  to  a  particular  study  period.  These  forecast  fields  or  ‘history 
fields’,  on  the  model’s  sigma  levels  and  at  the  model’s  horizontal  grid  points,  are  a  first 
guess  to  the  analysis  program  for  the  next  COAMPS  forecast.  The  assimilation  cycle 
continues  until  history  fields  are  produced  for  the  simulation  study  period.  Thus,  the 
mesoscale  character  of  the  flow  is  generated  and  maintained  by  implementing  a  data 
assimilation  update  cycle. 
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NOGAPS  Cold  Start 


Multivariate  optimum 
interpolation  analysis 


tau  (hr) - ► 


Figure  4;  Schematic  representation  of  the  COAMPS  data  assimilation  update  cycle.  The  procedure 
begins  with  a  12-hour  NOGAPS  cold  start.  In  this  example,  the  cold  start  is  initiated  36  hours  before 
the  initiation  of  the  study  period  at  tau=0  h.  The  cold  start  is  followed  by  two  12-hour  data  assimilation 
updates  using  COAMPS  forecast  fields.  A  third  data  assimilation  update  initiates  the  24-hour  forecast 
for  the  simulation  study  period.  Boundary  conditions  during  the  simulation  are  supplied  by  NOGAPS 
forecasts  every  12  hours. 


3.2.  COAMPS  Directory  Structures 

The  COAMPS  code  can  be  run  on  several  different  UNIX  platforms  including  the 
CRAY,  SGI,  DEC  ALPHA,  and  HP  machines.  The  process  for  implementing  COAMPS 
involves  first  obtaining  the  five  COAMPS  tar  files.  Requests  for  COAMPS  must  be 
submitted  to  Dr.  Richard  Hodur,  email:  hodur@nrlmry.navy.mil.  Also  forward  a  copy 
of  the  request  to  the  COAMPS  system  administrator  Sue  Chen,  email: 
chen@nrlmry.navy.mil.  The  five  tar  files  are  named  coamps#.tar,  templates.tar, 
database.tar,  documenttar,  and  utility.tar where  the  '#’  represents  the  version  release 
number. 

The  user  first  creates  a  /COAMPS  subdirectory  and  copies  the  COAMPS  tar  files 
into  it.  Then  each  of  the  tar  files  is  unarchived  using  the  UNIX  command:  “tar  -xvf  tar 
filenamd’.  Figure  5  shows  the  resulting  directory  structure  obtained  by  performing  this 
step,  along  with  several  additional  subdirectories  (marked  with  an  asterisk)  that  are 
automatically  created  during  the  process  of  setting  up  a  COAMPS  job. 

The  initial  COAMPS  directory  contains  five  subdirectories  comprising  the  entire 
COAMPS  system.  Referring  to  Figure  5,  these  subdirectories  are  listed  and  described 
below: 

COAMPS  Subdirectories _ 

•  /coamps#  -  COAMPS  master  source  codes,  libraries,  executables  and  prologues  for 
a  particular  version  release  number  given  here  by  the  '#’ 

/Makefile  -  UNIX  file  that  creates  (1)  library  archives  of  each  of  the  codes  in  /libsrc, 
and  (2)  an  executable  by  linking  the  libraries  with  the  driver  programs  in  /src 
/libsrc  -  source  codes  for  COAMPS  libraries 

/Makefile  -  UNIX  file  that  creates  a  library  archive  for  each  of  the  /libsrc  codes 
/coampslib  -  subdirectory  containing  COAMPS  source  code 
/fishpaklib  -  subdirectory  containing  NCAR  direct  solvers 
/fnoclib  -  subdirectory  containing  FNMOC  system  files 
/nl_beqlib  -  subdirectory  containing  nonlinear  balancing  code 
/oilib  -  subdirectory  containing  multivariate  optimal  interpolation  code 
/src  -  source  code  for  COAMPS  driver  programs 

/Makefile  -  UNIX  file  that  creates  the  COAMPS  executables 
/coamps_analysis  -  subdirectory  containing  the  analysis  driver  program 
/coamps_forecast  -  subdirectory  containing  the  forecast  driver  program 
/newdtg  -  subdirectory  containing  the  date-time  group  program 
/prologues  -  log  files  for  describing  source  codes  and  tracking  code  modifications 
^he  following  additional  subdirectories  {/lib,  /bin)  are  created  by  the 
Makefile  in  subdirectory  /coamps^ 

*/lib  -  location  of  library  archives  from  each  of  the  codes  in  /libsrc 
*/bin  -  location  of  executables  for  COAMPS  master  code 
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Figure  5:  CO  AMPS  directory  structure  contains  five  standard  subdirectories:  /coamps#,  \vhere  #  is  the  version 
release  number,  /database,  /templates,  /document,  and /utility.  These  subdirectories  are  obtained  when  the 
CO  AMPS  tar  files  are  unarchived.  The  two  additional  subdirectories,  /run  and  /mod  are  labeled 
with  an  asterisk  to  indicate  that  they  are  created  later  by  running  the  script  file  get.templates  in  the. /templates 
subdirectory.  Similarly,  the /coampsWlib  and /coamps#/bin  are  also  created  later  when  the  user  performs  a 
“make”  of  the  master  COAMPS  code  in  the  /coamps#  subdirectory. 


•  /templates  -  files  for  generating  templates  to  set  up  and  run  a  COAMPS  job 

/get.templates  -  script  file  for  creating  (1)  a  machine  dependent  run  script  used  to 
submit  a  COAMPS  job  and  (2)  a  standardized  COAMPS  makefile  used  to 
compile  modified  COAMPS  code 

•  /database  -  COAMPS  databases 

/masclim  -  subdirectory  containing  global  surface  climatology  database 
/masgiss  -  subdirectory  containing  Goddard  Institute  for  Space  Studies  database 
/masdted  -  subdirectory  containing  1  km  terrain  database 
/nogaps  -  subdirectory  containing  a  benchmark  set  of  NOGAPS 
/adp  -  subdirectory  containing  a  benchmaii<  set  of  observational  data 

•  /document  -  COAMPS  documentation 

/Users_guide  -  COAMPS  User’s  Guide 

•  /utility  -  COAMPS  graphics  and  utility  programs 

/Templates  -  script  files  for  processing  COAMPS  output 
(The  following  additional  subdirectories  {/mod/Templates,  /run/Templates)  are 
created  by  the  gettemplates  script  file) 

•  */run  -  run  scripts  for  and  output  from  COAMPS  case  studies 

*/Templates  -  standardized  run  scripts  for  submitting  COAMPS  jobs 

/run.coamps  -  run  script  template  created  by  running  get.templates 

•  */mod  -  modifications  to  COAMPS  source  code 

*/T emplates  -  standardized  makefile  for  creating  a  COAMPS  jobs 

/Makefile  -  makefile  template  created  by  running  get.templates 


Once  the  COAMPS  directory  structure  has  been  set  up,  the  code  must  be 
compiled.  This  step  requires  a  FORTRAN  77  or  FORTRAN  90  compiler  for  dynamic 
storage  allocation,  in  the  /coamps#  subdirectory,  type  “make  machine  namd’.  The 
‘machine  name’  can  be  chosen  from  the  following  options: 


Tcible  3  : 


Standard  machine  settings  for  COAMPS  compilation. 


MACHINE  NAME 

DESCRIPTION 

Cray 

CRAY 

crav_withjsis 

CRAY  using  operational  database  ISIS 

sgi 

SGI  R10000 

sgLg 

SGI  R10000  debug  mode 

alpha 

DEC  ALPHA 

alpha  q 

DEC  ALPHA  debug  mode 

hp 

HP  9000 

hp-g _ 

HP  9000  debug  mode 

This  command  creates  the  /coampsWIib  subdirectory  and  the  COAMPS  library  archives 
{coampsUb.a,  fishpacklib.a,  fnoclib.a,  nLbeqlib.a,  oilib.a)  as  well  as  the  /coamps#/bin 
subdirectory  and  the  master  code  executables  {coamps_analysis.exe, 
coamps_forecast.  exe) . 
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The  generation  of  these  COAMPS  executables  allows  the  user  to  proceed  to  the 
next  step:  running  a  COAMPS  job.  This  is  typically  accomplished  by  submitting  a 
COAMPS  run  script.  Obtaining  a  run  script  requires  running  the  script  file  get.templates 
located  in  the  /coamps#/templates  subdirectory.  Within  this  subdirectory,  type 
“get.templates  machine  name”  inserting  the  appropriate  machine  name  from  Table  3. 
This  command  prompts  the  user  for  input  and  upon  completion,  creates  two  new 
subdirectories;  subdirectory /coa/r7ps#/ru/7/re/r7p/afes  containing  a  standardized  run 
script  template  called  run.coamps,  and  subdirectory  /coamps#/mod/T amplates 
containing  a  standardized  makefile  template  called  Makefile.  This  Makefile  is  used  to 
create  updated  COAMPS  executables  when  the  user  modifies  the  original  master 
source  code.  The  procedure  for  modifying  the  master  COAMPS  code  and  submitting  a 
COAMPS  simulation  is  discussed  in  Section  4.  First,  a  review  of  the  COAMPS  coding 
framework  is  given  in  the  next  section, 

3.3.  COAMPS  Code 

COAMPS  consists  of  two  separate  driver  programs  that  execute  the  analysis  and 
the  forecast.  In  the  following  sections,  the  flow  chart  diagrams  and  primary  subroutines 
associated  with  each  of  these  two  programs  are  summarized  and  outlined. 

3.3.1 .  Analysis 

The  purpose  of  the  analysis  routine  coama.F  is  to  produce  a  set  of  initial 
conditions  for  the  forecast  model.  In  addition  to  specifying  the  surface  parameters  and 
basic  state  arrays,  the  analysis  specifies  3D  atmospheric  fields  either  by  interpolating 
NOGAPS  fields  to  the  COAMPS  grids  or  by  blending  observations  with  the  first  guess 
fields  using  the  MVOI  (multivariate  optimum  interpolation)  analysis.  Additionally, 
coama.F  prepares  lateral  boundary  conditions  for  the  forecast  model. 

Figure  6  shows  a  diagram  of  the  analysis  routine’s  input/output  including  user- 
defined  namelist  values,  first  guess  fields,  and  surface  databases  as  input,  and  initial 
forecast  fields  as  output.  Each  of  the  inputs  are  specified  by  the  user  to  customize  the 
initial  conditions  and  properly  set  up  a  COAMPS  simulation.  In  the  analysis  code,  the 
gridnl  namelist  parameters  are  passed  in  as  subroutine  arguments.  These  parameters 
include:  (1 )  the  grid  dimensioning  variables,  which  are  used  later  in  the  code  to  indicate 
do-loop  limits  (ma(nn),na(nn),kka),  and  (2)  the  array  allocation  length  nwords  which  is 
used  to  dimension  the  total  analysis  program  array  a(nwords).  Additionally,  the  pointer 
location  arrays  (ivar  (nn))  are  dimensioned  in  common  blocks  which  are  incorporated 
via  the  include  file  apointers.h.  Then  the  coamni  namelist  parameters  are  read  In  to 
specify  model  features  and  case  study  details.  This  namelist,  containing  nearly  200 
input  variables,  is  the  primary  vehicle  for  the  user  to  easily  interact  with  the  system 
without  actually  making  coding  changes.  The  following  information  is  provided  to  the 
analysis  routine  through  the  coamn/ namelist  inputs: 
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INPUT 


coama.f  I/O 


OUTPUT 


I 


USER  INPUT 


[-  gridnlmmelist  - 

I-  coamnl  namelist  - 

•  INITIAL  DATA 

h  •  SURFACE  FIELDS 

1-  observational  databases 


(NAMELIST  DIRS:  dsclim,  dsgiss,  dsdted) 
I—  user  defined  idealized  surface  - 


E 


real  data 


(NAMELIST  VARS;  seatmp,  albdo,  zrough,  sfcwet,  alndpct) 

•  ATMOSPHERIC  FIELDS 

Yit —  user  defined  idealized  data - ►I 

r  jcajjQ-s^NAMELIST  ARYS;  psnd,  qsnd,  tsnd,  usnd,  vsnd) 

>  Oy 

Operational 
“ISSdatabase 

^  (NAMELIST  DIR:  dsetng) 

fields . 

R&D 

■flat  files . . 

(NAMELIST  DIR:  dsngfl) 

Operational 
IMS  data  base  ■ 


COAMPS 

fields 


observational  data 


■(NAMELIST  DIRS:  dsetg) 


R&D 
flat  files  ' 


(NAMELIST  DIR:  dsnrft) 


(NAMELIST  DIR:  datfil) 


STANDARD  OUTPUT/ERROR  FILE  - ► 


IDEALIZED  FIELDS 

(NAMm.ISTDIR:dsnrff) 


REAL  DATA 
ANALYZED 
FIELDS 


C^rational 
IMS  data  base- 


(NAMELISTDIR:  dsetg) 


R&D 

flat  files  — — .  ► 

(NAMELIST  DIR:  dsnrff) 


Figure  6:  Flow  chart  of  input/output  for  the  analysis  routine  coama.f.  The  analysis  routine’s  input  requirements 
include  user  input  from  the  namelists,  as  well  as  initial  surface  and  atmospheric  fields.  The  analysis  routine  outputs 
a  standard  output/error  file  and  the  idealized  or  analyzed  real  data  fields  for  the  forecast  program.  For  reference, 
corresponding  namelist  parameters  are  shown  in  parentheses. 


Type  of 

Model  Run  (Namelist  variables) 

Directory  Locations 

of  initial  Input  Data  (Namelist  variables) 

•operational  or  R&D  (Idbms,  Idbmi,  Idbmo) 
•idealized  or  real  data  (icase) 

•data  assimilation  update  (iupd,  loi,  loimf) 

•surface  databases  (dsclim,  dsgiss,  dsdted) 
•observational  data  (datfil) 

•NOGAPS  fields  (dsetng,  dsngff) 

•COAMPS  fields  (dsetg,  dsnrff) 

These  parameters,  defined  further  in  Appendix  D,  indicate  the  location  and  type  of  data 
required  to  set  up  the  initial  conditions.  Note  that  before  the  analysis  can  be  run,  the 
user  must  obtain,  and  specify  the  directory  location  of  this  initial  input  data.  Either 
idealized  or  real  data  may  be  used  to  provide  the  initial  fields.  Furthermore,  real  data 
model  runs  can  be  done  in  either  an  operational  or  R&D  mode. 

Operational  use  of  COAMPS  is  accomplished  by  toggling  the  three  database 
logical  variables  (idbms,  Idbmi,  Idbmo)  to  true,  accessing  Fleet  Numerical 
Meteorology  and  Oceanography  Centers  (FNMOC)  operational  databases  that  are 
managed  by  the  Integrated  Stored  Information  System  (ISIS).  For  R&D  purposes 
however,  COAMPS  reads  in  real  data  initial  input  fields  from  ‘flat  files’  using  32-bit  IEEE 
unformatted  FORTRAN  I/O.  In  that  case,  the  user  must  obtain  the  required  data  (first 
guess  fields,  surface  databases,  observations  and  lateral  boundary  conditions),  in  the 
required  format  (32-bit  IEEE  unformatted  FORTRAN  files),  and  specify  the  directory 
location  in  which  each  resides  (dsclim,  dsgiss,  dsdted,  datfil,  dsetng,  dsngff,  dsetg, 
dsnrff).  Once  the  initial  input  fields  are  processed  in  coama.F,  they  are  written  out  to 
either  ‘restart’ files  for  idealized  R&D  output,  to  the  ISIS  database  for  operational 
output,  or  to  flat  files  for  real  data  R&D  output.  The  flat  files  are  named  according  to  a 
standardized  naming  convention  given  in  Section  3.3.2.3  Writing  Output,  Table  4. 

Proper  processing  of  the  initial  input  fields  to  obtain  initial  conditions  fields 
requires  an  understanding  of  the  coama.F  coding  structure.  The  flow  diagram  shown  in 
Figure  7  lists  the  subroutines  called  by  coama.F.  The  first  section  of  the  analysis 
routine  consists  of  reading  in  the  namelist  and  performing  initial  setup  functions.  Once 
the  coamni and  gr/dr?/ namelist  are  read  in,  calls  to  subroutines  chekmn.F,  checkz.F, 
grid.F  and  rtm^t/v.F specify  the  model  domain  space,  map  factors  and  grid  nesting. 

Then  the  surface  characteristics  for  each  nest  are  determined  from  subroutine 
calls  to  sfcpar.F,  getsst.F,  user_sfc.F,  and  soUtp.F .  In  sfcpar.F,  the  surface  fields 
(riand,  albedo,  gwet,  zO,.csst)  are  initially  read  in  from  a  global  climatology  database. 
For  specific  regions,  higher  resolution  data  for  the  albedo  and  zO,  available  from  the 
Goddard  Institute  for  Space  Studies  (GISS),  replaces  the  climatological  values.  In 
subroutine  getsst.F,  the  resolution  of  the  surface  temperature  tsea  over  water  is 
improved  by  using  a  1°  NOGAPS  field  or  a  125km  x  125km  hemispheric  grid  derived 
from  the  Optimum  Thermal  Interpolation  System  (OTIS),  or  from  COAMPS  history 
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setup  grid 

chekmn.f,  checkz.f,  grid.f,  hm2uv.f: 
setup  model  domain,  map  factors,  grid  nesting 


user  input 

coamni 

namelist 


setup  surface  and  terrain  fields 

sfcpar.f,  getsstf,  user_sfc.f,  soiltp.f, 
atopo.f,  tmatch.f 

assign  surface  and  terrain  fields  from  databases  or 
from  horizontally  homogeneous  user  defined  values 


read  and  process  initial  input  fields 


icase^iO. 


freturi 


idealized  data 

specify  fields 

refsndJ,  istate.f,  flds§#.f 
create  mean  state  and  initial  idealized 
fields  from  reference  sounding;  add 
perturbafons  to  initial  state 


interpolate  fields 

icm2fg.f 

interpolate  coarse  mesh  (1)  fields  to 
inner  nests  for  inhomogeneous  initial 
conditions 


_ write  fields _ 

idealized  data 

ioanl.f  (iozavg.f,  iosfcO.f, 
iosfctfjoatmt.f) 
write  idealized  initial  fields  to 
‘resfarf’ files 

real  data 

outhty.f,  iosfc.f,  outanl.f,  outinc.f 
write  surface  fields,  analyzed  fields 
and  MVOI  increments 


real  data 

specify  mean  state 

astate.f 

create  mean  state  from  standard  atmosphere 

update  surface  fields 

gethty.f,  nethty.f 

update  surface  fields  from  COAMPS  history 
fields 

get  first  guess  fields 

getngJ,  ngSfg.f,  fg2fg.f,  getfgl.f 
read  first  guess  fields  on  pressure  levels  from 
(1)  NOGAPS  (interpolate  to  COAMPS  grids),  or 
(2)  COAMPS  history  fields 

perform  multivariate  optimum 
interpolation  analysis 

oi_anl.f,  anlfldJ,  netinc.f 
incorporate  observations 

get  boundary  conditions 

getbdy.f,  tendbd.f,  getng.f 
extract  coarse  mesh  boundary  values  from 
NOQAPS  fields  for  the  remainder  of  the  study 
period  and  write  out  tendancies 


Figure  7 :  Flow  chart  of  the  analysis  routine  coama.f.  The  analysis  routine  is  described  in  five 
parts:  setup  parameters,  setup  grid,  setup  surface  and  terrain  fields,  read  and  process  initial  input 
fields  and  write  fields.  The  subroutines  called  in  each  part  are  shown  in  italics.  The  atmospheric 
fields  may  be  either  idealized  or  real  data  as  indicated  by  the  value  of  namelist  parameter  lease. 
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field  if  either  of  these  source  are  available.  Otherwise,  the  climagologicai  value  csst  is 
used,  in  subroutine  soiltp.F,  values  for  the  deep  soil  temperature  array  tsoil  are 
generated  from  climatological  data.  Also  generated  in  subroutine  sfcpar.F  is  the 
surface  topography  field  (zsfc),  initialized  from  a  20km  terrain  dataset.  In  certain 
regions  around  the  globe,  1  km  resolution  terrain  data  is  used  for  meshes  with  grid 
spacing  less  than  20km. 

In  setting  up  the  surface  characteristics  of  a  particular  case  study,  subroutine 
i/se/isfc.F allows  the  user  to  customize  the  surface  parameters  with  idealized  values 
through  namelist  input.  For  example,  user-defined,  horizontally  homogeneous  values 
for  the  land/sea  table,  surface  temperature,  roughness,  albedo  and  ground  wetness 
arrays  are  given  by  namelist  parameters  aindpct,  seatmp,  zrough,  albdo,  sfcwet 
when  the  namelist  flags  iindfig,  iseafig,  izOfig,  iaibfig,  igwtfig  are  set  to  ‘O’.  The 
namelist  flag  itopofig  when  set  to  ‘0’  assigns  the  terrain  height  to  zero  everywhere. 
Finally,  subroutine  fmafch.F  matches  the  terrain  field  across  the  mesh  boundaries  and 
subroutine  afopo.F computes  the  terrain  gradient  arrays.  The  coding  then  branches  for 
idealized  or  real  data  initial  conditions  as  described  in  the  next  two  sections. 

3.3.1 .1.  Idealized  Data 

Idealized  initial  conditions  are  typically  provided  by  user-defined  profiles  of  winds, 
temperature,  pressure  and  moisture  assumed  to  be  horizontally  homogeneous  across 
the  model  domain.  These  profiles  are  read  in  through  namelist  arrays  usnd,  vsnd, 
tsnd,  psnd,  and  qsnd  which  are  used  when  the  namelist  parameter  lease  is  nonzero. 
The  subroutine  refenc/.F  converts  the  initial  user-defined  reference  sounding,  into  MKS 
units.  Then  a  homogeneous  idealized  mean  state  is  determined  from  that  data  in 
subroutine  istate.F  for  each  of  the  model  nests.  Perturbations  can  be  added  to  the 
mean  state  by  specifying  a  value  for  the  namelist  parameter  icase=’##’.  This  number 
'##’,  corresponds  to  a  subroutine  f/afs##.F  where  the  user  assigns  the  field 
perturbations.  When  using  inhomogeneous  idealized  fields,  set  the  namelist  parameter 
jcm2fg  to  1  to  interpolate  the  outer  coarse  mesh  (1 )  fields  to  the  inner  nests.  This 
interpolation  is  done  in  subroutine  icm2fg.F.  In  subroutine  ioanl.F,  the  complete  set  of 
idealized  Initial  conditions  (surface,  basic  state  and  prognostic  fields)  is  written  to  files 
that  begin  with  the  prefix  ‘restart.  These  ‘restart  files  become  the  initial  forecast  fields 
read  in  by  coamm.FXo  begin  an  idealized  simulation.  As  discussed  later,  the  ‘restart 
files  can  be  created  at  any  time  during  the  COAMPS  model  run  and  are  also  used  to 
initiate  COAMPS  at  a  nonzero  forecast  time  (ktaust^tO). 

3.3.1 .2.  Real  Data 

An  alternative  to  idealized  initial  input  fields  is  the  use  of  global  or  mesoscale  first 
guess  fields  that  may  be  updated  with  observations.  The  first  guess  can  be  obtained 
from  either  a  NOGAPS  analysis  or  forecast,  or  from  a  previous  COAMPS  forecast  valid 
at  the  desired  date-time  group  (history  fields). 
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First,  the  basic  state  profiles  are  specified  through  subroutine  astate.F .  Then 
the  previously  specified  surface  parameters  are  replaced  with  higher  resolution 
COAMPS  history  fields,  if  they  are  available.  From  calls  to  gethty.Fan6  nethty.F,  the 
2D  arrays  for  surface  parameters  on  each  nest  are  updated  if  the  logical  namelist 
parameter  lnrhty=.true..  These  fields  (gwet,  zO,  tsea),  Initially  read  in  from  a  global 
climatology  database,  are  assigned  values  given  by  the  COAMPS  history  fields.  While 
the  array  for  snow  coverage  (snow)  is  given  values  from  the  NOGAPS  global  field,  the 
other  2D  surface  fields  obtain  values  from  the  COAMPS  history  fields.  These  fields 
include:  the  10m  winds  (ulOm,  vIOm),  boundary  layer  depth  (biht),  surface  fluxes 
(hfixs,  hfixl)  and  wind  stress  (stres).  When  lnrhty=.false.,  or  if  history  fields  are  not 
found,  this  latter  set  of  fields  are  set  to  zero  initially. 

The  procedure  for  obtaining  the  first  guess  fields  begins  with  subroutine  getng.F 
that  reads  in  NOGAPS  fields  for: 

•  surface  pressure 

•  u-  and  v-momentum 

•  geopotential  heights 

•  vapor  pressure 

These  fields  are  specified  at  1°  horizontal  resolution  and  at  pressure  levels  pr(lm), 
where  Im  and  pr  are  user  defined  namelist  parameters  defined  further  in  Appendix  D. 
The  fields  are  then  interpolated  to  the  coarse  mesh  (1)  grid  points  in  subroutine 
ng2fg.F.  Subroutine  fgr2fgr.F performs  the  interpolation  of  the  coarse  mesh  (1)  to  the 
inner  nests.  If  a  previous  COAMPS  forecast  is  available  and  the  proper  namelist  flag  is 
set  (Iupd?i0),  then  subroutine  getfg1.F  \s  called  to  overwrite  the  above  fields  with 
COAMPS  history  fields  for  each  nest.  Both  routines  getng.F  and  getfgl.F  search 
backward  through  seven  date-time  groups  in  an  attempt  to  provide  COAMPS  with  the 
most  recent  data. 

Finally,  the  first  guess  fields  are  adjusted  based  upon  observational  data  via  a 
multivariate  optimum  interpolation  (MVOI)  analysis  performed  in  the  subroutine  oi_anl.F 
(lol=.true.).  Note  that  the  MVOI  analysis  subroutines  are  not  stored  with  the  rest  of  the 
COAMPS  source  code  (See  Section  3.2,  and  Figure  5),  but  are  contained  in  a  separate 
subdirectory  and  compiled  into  a  separate  library.  The  o/Lan/.F subroutine  determines 
height  and  momentum  increments  based  upon  differences  between  the  first  guess 
fields  on  the  coarse  mesh  (1)  and  ADP  observations  (i.e.  rawinsonde,  SSMI,  satellite 
derived  data,  etc.).  These  increments  are  added  to  the  first  guess  fields  and  a  9-point 
smoother  is  applied  in  subroutine  anIfId.F.  The  same  procedure  can  be  used  to  obtain 
MVOI  increments  directly  on  the  inner  nest  first  guess  fields  (loimf=.true.).  Alternatively 
(loimf=.false.),  the  coarse  mesh  (1)  MVOI  increments  can  be  interpolated  to  the  inner 
nest  grids  and  then  added  to  the  inner  nest  first  guess  fields  using  subroutine  netinc.F] 
however,  this  procedure  is  less  effective  at  maintaining  the  mesoscale  structure  of  the 
observations.  Using  a  Cressman  scheme,  additional  analyses  are  also  available  on  the 


24 


pressure  levels  for  temperature  (ltanl=.true.)  and  dew  point  depression  (lqanl=.true.), 
as  well  as  for  surface  temperature  (ltanls=.true.). 

To  produce  a  forecast  from  real  data,  boundary  conditions  for  the  prognostic 
fields  must  be  supplied  to  the  coarse  mesh  (1)  at  a  given  time  interval  (itauin).  For  a 
specified  forecast  length  (itauf),  the  boundary  values  are  extracted  from  the  NOGAPS 
forecast  fields  in  subroutine  getbdy.F,  and  the  tendencies  are  computed  and  written  out 
in  subroutine  tendbd.F,  again  following  the  standard  naming  convention. 

The  final  set  of  real  data  initial  conditions  is  written  out  by  three  output 
subroutines.  Subroutines  outhty.F,  /osfc.Fwrite  out  the  surface  fields  while  subroutine 
oufan/.F  writes  out  the  analyzed  fields  at  the  pr(lm)  pressure  levels.  For  use  in  data 
assimilation,  subroutine  outinc.F  writes  out  the  MVOI  increments  (loi=.true.).  Since  the 
basic  state  arrays  are  easily  recomputed,  they  are  not  written  out  by  the  analysis 
routine  for  real  data  cases.  The  files  created  by  the  analysis  are  named  following  the 
standard  convention  described  in  Section  3.3.2.3  for  which  each  filename  begins  with  a 
four  character  prefix  identifying  the  field  of  data  it  contains.  This  completes  the  coding 
structure  associated  with  the  analysis  subroutine  coama.F and  the  preparation  of 
COAMPS  analyzed  fields. 

3.3.2.  Forecast 

The  purpose  of  the  coamm.F  code  is  to  produce  predicted  values  of  the  time- 
dependent  variables.  Associated  input/output  tasks  include  reading  in  the  initial 
conditions  and  writing  out  the  model  results.  These  I/O  functions  are  shown  in  Figure  8. 
As  in  the  analysis  routine,  first  the  gr/c/n/ namelist  parameters,  passed  in  as  subroutine 
arguments,  are  used  to  indicate  do-loop  limits  (ma(nn),na(nn),kka),  and  the  array 
allocation  length  nwords  is  used  to  dimension  the  total  forecast  program  array 
a(nwords).  Additionally,  the  pointer  location  arrays  are  dimensioned  in  common  blocks 
which  are  incorporated  via  the  include  file  mpointers.h.  Then  the  coam/?/ namelist 
parameters  are  read  in  to  indicate  model  features  and  case  study  details  specified  by 
the  user.  The  2D  graphics  instructions,  designated  here  by  filenames  ‘ocards’  and 
‘xcardd,  are  read  in  so  that  the  proper  output  fields  will  be  saved  for  graphical  display. 
Finally,  the  analysis  fields  created  by  coama.F  are  input  from  either  the  ‘restart' files, 
ISIS  database,  or  R&D  flat  files.  The  forecast  routine  produces  output  at  specified  time 
intervals  for  evaluation  of  model  results.  COAMPS  output  includes  a  standard 
output/error  file,  and  additional  output  as  requested  by  the  user:  ‘restart  files,  1 D,  2D 
and  3D  graphics  files. 

With  an  understanding  of  the  forecast  routines  I/O,  the  coamm.F  coding  is 
described  below.  The  forecast  is  produced  by  integrating  a  set  of  model  equations 
governing  the  prognostic  nature  of  each  time-dependent  variable.  Where  necessary,  to 
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Figure  8;  Flow  chart  of  input/output  for  the  forecast  routine  coamm.f.  The  forecast  routine  uses  from  the 
coamnl  namelist,  input  parameters  lease  to  setup  idealized  or  real  data,  ktaust  to  determine  if  'restart'  files 
or  if  the  analysis’  first  guess  fields  are  input,  and  idbms  to  access  the  operational  database  or  R  &  D 
flat  files.  The  coamnl  namelist  also  assigns  the  directory  locations,  shown  in  parentheses,  of  data  input/output. 
The  files  ocards  and  xcards  contain  user-specified  graphics  directives.  For  reference,  corresponding  namelist 
parameters  are  shown  in  parentheses. 


close  the  system  of  equations  or  to  represent  certain  processes,  physical 
parameterizations  have  been  utilized.  The  following  section  discusses  the  coding 
structure  and  subroutine  calls  within  subroutine  coamm.F.  Aside  from  initial  setup 
functions  performed  by  the  subroutine,  it  is  subdivided  into  two  other  parts:  reading  in 
initial  conditions  and  performing  the  forecast.  While  reviewing  this  section,  refer  to  the 
flow  chart  and  subroutines  shown  n  Figure  9. 

3.3.2.1.  Reading  Initial  Conditions 

Subroutine  coamm.F  begins  by  reading  in  the  coamn/ namelist  to  allow  for  user- 
defined  model  input  and  then  performs  several  setup  functions  that  initialize  and  define 
constants,  parameters,  variables  and  arrays  used  later  in  the  forecast  routine. 

Additional  setup  code  is  embedded  in  subroutine  coami.F:  the  model  domain  space, 
map  factors  and  grid  nesting  are  determined  from  routines  grid.F,  hm2uv.F,  and  2D 
graphics  instructions  are  read  in  from  routines  reado.F,  readx.F.  Then,  the  assignment 
of  COAMPS  initial  forecast  fields  is  done  in  coami.F.  If  idealized  data  (icase?K))  are 
used  or  a  restart  run  is  initiated  (ktaust^^O),  subroutine  iomdl.Fls  called  to  open  the 
'restart  files  to  provide  the  initial  conditions.  Otherwise,  real  data  initial  conditions  are 
obtained  from  subroutine  calls  to  insfc.F,  atopo.F,  astate.F,  inIvl.F. 

When  using  real  data,  subroutine  /nMF  specifies  the  initial  forecast  fields  based 
upon  the  type  of  data  assimilation  select  by  the  namelist  parameter  iupd.  When 
iupd=0,  a  NOGAPS  cold  start  is  performed  by  interpolating  analyzed  NOGAPS  fields 
on  pressure  levels  to  COAMPS  sigma  levels  and  horizontal  grid  points.  When  iupd=1 , 
a  full  COAMPS  update  is  performed  by  interpolating  the  analyzed  COAMPS  fields  on 
pressure  levels  to  sigma  levels  (subroutines  instdp.F,  stdp2z.F).  When  iupd=2,  an 
incremental  update  is  performed  in  subroutine  incrup.F,  where  the  COAMPS  history 
fields  on  model  sigma  levels  and  the  MVOI  increments  on  pressure  levels  are  read  in. 
After  the  increments  are  interpolated  to  the  model  sigma  levels,  variational  adjustments 
are  made  on  the  pressure  and  potential  temperature  increments  to  bring  the  fields  into 
hydrostatic  balance  and  then  the  increments  are  added  to  the  COAMPS  forecast. 

Returning  back  to  subroutine  coami.F,  next  a  few  additional  preliminary  steps  are 
taken  before  performing  the  model  forecast.  Boundary  values  are  extracted  in 
subroutines  readbd.F an6  tendbd.F,  and  several  model  parameters  are  computed  using 
namelist  input.  For  example,  lateral  boundaries  and  weighting  functions  are  determined 
from  kgetbc,  ibdya,  Jbdya;  the  mixing  length  array  is  defined  using  iamxgl;  the  sponge 
layer  is  specified  by  Iralee,  Ispong,  nrdamp;  and  the  diffusion  coefficients  are 
determined  from  dif4th  and  dif2nd.  In  addition,  the  prognostic  fields  are  broadcast  to 
the  other  time  levels  in  asetup.F.  If  the  logical  variable  linit=.true.,  subroutine  ainit.F  ls 
used  to  modify  or  update  certain  fields  used  by  the  model.  (Each  of  these  user- 
specified  parameters  are  further  defined  in  Appendix  D.)  Finally,  the  initial  forecast 
fields  are  written  out  by  subroutine  output. F. 
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stdp2z.f:  interpolate  pressure  level  fields 
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quick  print  data  to  standard  output,  1 D,  2D  and 
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Figure  9;  How  chart  of  the  forecast  routine  coamm.f.  The  forecast  routine  is  described  in 
five  parts:  setup  parameters,  setup  grid,  setup  2D  output,  read  initial  conditions  and  perform 
forecast.  The  subroutines  called  in  each  part  are  shown  in  italics.  Depending  on  the  value  of 
two  input  parameters,  lease  and  ktaust,  co<3mm./reads  data  from  the  ^restart’  files,  or  reads 
real  data.  For  real  data  assimilation,  input  parameter  iupd  indicates  a  ‘full’  or  ‘incremental’ 
update.  Input  parameters  ktaust,  ktauf,  delta  determine  the  number  of  coarse  mesh 
iterations  (iterl).  Within  this  loop,  the  coarse  mesh  and  inner  nest  fields  are  integrated  in 
subroutine  amodel.f  (see  Figure  10).  Fields  necessary  to  restart  CO  AMPS  are  saved  when 
indicated  by  ksavea,  and  various  output  parameters  indicate  when  forecasted  data  is  written 
out  (see  Figure  11). 
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At  this  point,  control  returns  to  the  main  forecast  subroutine  coamm.F.  And  since 
all  model  arrays,  variables  and  parameters  have  been  assigned  values,  COAMPS  is 
now  ready  to  begin  integration  of  the  initial  forecast  fields. 

3.3.2.2.  Model  Integration:  Forecast  Loop  and  Atmospheric  Model 

Predicted  values  of  the  model  fields  are  obtained  by  second  order  leap  frog 
integration  of  the  nonhydrostatic  system  of  equations.  These  atmospheric  equations 
are  coded  in  subroutine  amode/.F  which  is  discussed  in  the  next  paragraph.  The 
integration  of  the  prognostic  fields  is  controlled  by  ‘the  forecast  loop’  of  coamm.F  (See 
Fig.  9).  The  loop  indices  iterls  and  iterle  are  determined  from  user-specified  namelist 
input  ktaua  and  ktauf.  Each  pass  through  the  forecast  loop  integrates  the  coarse 
mesh  (1)  fields  forward  one  delta  time  step.  In  real  data  cases,  coarse  mesh  (1) 
boundary  values  are  fTrst  obtained  from  subroutine  readbd.F every  kgetbc  hours.  Then 
the  predicted  (time  level  3)  coarse  mesh  (1)  fields  are  returned  by  subroutine  amodel.F. 
If  parameter  nnest>1 ,  this  step  is  followed  by  nested  calls  to  amodel.Fto  integrate  the 
fields  on  each  of  the  Inner  nests.  Because  each  inner  nest  is  required  to  be  one-third 
the  horizontal  grid  spacing  of  it’s  parent  mesh,  it  is  iterated  three  times  for  each  parent 
mesh  iteration,  over  a  time  step  that  is  one-third  that  of  it’s  parent .  At  the  beginning  of 
each  inner  nest  iteration  loop,  boundary  values  are  extracted  from  the  parent  mesh  in 
subroutine  bdr.F and  mbndy.F  (for  real  data  cases  and  only  when  nbdya=7)  and  the 
predicted  inner  nest  fields  are  returned  by  subroutine  amodel.F.  At  the  end  of  each 
inner  nest  iteration  loop,  if  namelist  parameter  I2way=.true.,  the  forecasted  fields  are 
fed  back  to  the  parent  mesh  in  subroutine  feedbk.FXo  allow  for  a  two-way  interactive 
grids.  Before  continuing  to  the  next  coarse  mesh  (1 )  time  step,  the  forecast  time  is 
checked  to  see  if  it  coincides  with  times  specified  by  the  user  for  saving  ‘resfa/t' files 
(ksavea)  or  model  output  (see  Figure  1 1  and  Section  3.3.2.3).  Of  primary  importance 
in  the  forecast  loop  is  the  call  to  the  atmospheric  model  subroutine  amode/. F  described 
below. 


In  general,  all  namelist  inputs,  model  parameters  or  arrays  used  during  the 
model  integration  are  passed  through  the  calling  arguments  to  amode/.F  rather  than 
through  common  blocks.  This  allows  for  more  shared  memory  usage  and  dynamic 
memory  allocation  within  each  subroutine.  As  the  subroutines  called  by  amodel.F  are 
described,  refer  to  the  atmospheric  model  flow  chart  in  Figure  10.  The  subroutine  calls 
within  amodel.F  are  separated  into  four  categories: 

•  diagnostics  -  asigdl.F,  adIvSd.F,  cgterm.F,  adiagk.F 

•  parameterizations  -  radiat.F,  sfcflx.F,  afore.F,  kfpara.F,  adjtq.F 

•  prognostic  equations  -  afore.F,  aforth.F,  arhsu.F,  arhsv.F,  arhsw.F,  arhsp.F,  alhs.F 

•  adjustments  -  ashift.F,  abdye.F,  abdytq.F,  abdyuv.F,  abdyp.F,  amixtq.F,  amixnf.F, 

amixuv.F,  amixw.F,  atimef.F 
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amodeif 

atmospheric  model  equations 
and  parameterizations 


ashift.f  (adj) :  shift  time  level  arrays  to  prepare  for  the  next  iteration 
{e1,  thi,  qvl,  qcl,  qrl,  qil,  qsl,  u1,  v1,  w1,  p1, 
e2,  th2,  qv2,  qc2,  qr2,  qi2,  qs2,  u2,  v2,  w2,  p2, 
e3,  th3,  qv3,  qc3,  qr3,  qi3,  qs3,  u3,  v3,  w3,  p3) 


asigdtf  (diag):  compute  the  vertical  velocity  associated  with  the 
model  coordinates  for  time  level  2  (sigdt) 


adiv3d.f (diag):  compute 
the  3D  divergence  array 
(div3) 


cgterm.  f  (diag) :  compute 
the  counter  gradient  terms 
(th3,  qv3) 


adiagkJ  (dAag):  compute 
the  Smagorinsky  eddy  mixing 
coefficient  array  (eks) 


sfcfix.  f  (para) :  compute  the 
surface  fluxes  and  stresses 
(hfluxi,  hfluxs,  stresx,  stresy) 


afore.f  (para/prog):  predict  tendencies  and  update  the  turbulent  kinetic  energy  (e3) 


aforth.f  (prog):  predict  tendencies  for  the  potential  temperature  array  (thrhs) 


aforqx.f  (prog):  predict 

tendencies  of  cloud  and  rain 

mixing  ratios  (qc3,  qr3) 

aforqx.f  (prog):  predict  tendencies  of  ice  and  snow  mixing  ratios  (qi3,  qs3) 


aforqx.f  (prog):  predict  tendencies  of  water  vapor  mixing  ratio  (qv3) 


arhsu.f  {prog):  predict  advective  mode  tendencies  for  the  u-momentum  (urhs) 


arhsv.f  (prog):  predict  advective  mode  tendencies  for  the  v-momentum  (vrhs) 


arhsw.f  (prog):  predict  advective  mode  tendencies  for  the  vertical  Velocity  (wrhs) 


arhsp.f  (prog):  predict  advective  mode  tendencies  for  the  pressure  (prhs) 


abdye.f,  abdytq.f  (adj):  apply  boundary  conditions  (e3,  qc3,  qi3,  qr3,  qs3,  qv3,  th3) 


asigdtf  (diag):  compute  the  vertical  velocity  associated  with  the 
model  coordinates  for  time  level  2  (sigdt) 


(continued) 
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(continued) 


iter2  =  1,  ktaua:  mtaua  - - 


alhs.f  (prog/adj):  predict  the  fast  mode  tendencies  for  the  pressure, 
vertical  velocity  and  momentum  terms,  update  the  u-momentum  (u3), 
v-momentum  {v3),  vertical  velocity  (w3),  pressure  {p3),and  potential 
temperature  (th3),  apply  boundary  conditions  (th3,  p3,  u3,  v3,  w3) 


Irad  >EL 


update  th3  with  the  radiational  heating  rate  tendencies  (trad) 


Itke 
\  ?  y 


abdye.f,  abdytq.f,  abdyuv.f  (adj):  apply  boundary  conditions 
(e3,  th3,  qv3,  qc3,  qr3,  qi3,  qs3,  u3,  v3,  w3 


amixtq.f  (para):  implicit  vertical  mixing  (th3,  qv3) 


amixnf.f  (para):  implicit  vertical  mixing  (e3) 


amixuv.f  (para);  implicit  vertical  mixing  (u3,  v3) 


amixw.f  (para);  implicit  vertical  mixing  (w3) 


amixnf.f  {para):  implicit 
vertical  mixing  (qc3,  qr3) 


amixnf.f  (para):  implicit  vertical  mixing  (qi3,  qs3) 


(icupa^^ 

F\^'V3elxx\T 


abdye.f,  abdytq.f,  abdyuv.f  (adj):  apply  boundary  conditions 
(e3,  th3,  qv3,  qc3,  qr3,  qi3,  qs3,  u3,  v3,  w3 


kfpara.f  (para):  compute  convective  tendencies  {qc3,  qr3,  qv3,  th3) 


adjtq.f  (para):  compute  the  microphysical  tendencies  (th3,  qv3,  qc3,  qi3,  qr3,  qs3) 


radiat.f  (para);  compute  radiative  heating  rate  tendencies  (trad) 


abdytq.f  (adj):  apply  boundary  conditions  (th3,  qv3,  qc3,  qi3,  qr3,  qs3) 


atimef.f  (adj);  filter  high-frequency  temporal  oscillations 
(e2,  p2,  qc2,  qi2,  qr2,  qs2,  qv2,  th2,  u2,  v2,  w2) 


Figure  10:  Flow  chart  of  the  atmospheric  model  subroutine  Subroutines  called 

by  amodel.faxe  designated  as:  adjustment  (adj),  diagnostic  (diag),  prognostic  (prog),  and/or 
parameterization  (para)  routines.  Shown  in  parentheses  are  the  arrays  updated  in  each 
subroutinecalled  by  amodelf.  The  number  associated  with  each  array  indicates  the  time 
level.  Many  of  the  model  parameterization  features  are  controlled  by  input  parameters.  For 
example,  subgrid  scale  mixing  is  controlled  by  parameters  Itke,  iamxgl,  iahsgm,  surface 
fluxes  by  Iflux,  thermodynamic  processes  by  Imoist,  lice,  convective  parameterization  by 
Icupar,  dxmeso,  and  radiative  affects  by  Irad  and  the  internal  logical  parameter  luprad 
which  obtains  a  value  based  upon  the  values  of  both  Irad  and  dtrad.  See  Appendix  D  for  a 
more  thorough  description  of  these,  and  the  other  CO  AMPS  namelist  input  parameters. 
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Within  amodel.F,  first  the  time  levels  associated  with  the  prognostic  arrays  are 
shifted  in  subroutine  ashift.FXo  prepare  them  for  the  next  model  iteration.  To 
illustrate  this  process,  we  use  the  u-momentum  component  array  as  an  example: 

u1(m,n,kk)  =  u2(m,n,kk) 
u2(m,n,kk)  =  u3(m,n,kk) 
u3(m,n,kk)  =  u1(m,n,kk). 

After  ashift.F,  time  level  1  arrays  contain  previous  values  (time=t-At),  time  level  2  arrays 
contain  present  values  (time=t),  and  by  the  end  of  amodel.F,  time  level  3  arrays  contain 
predicted  values  (time=t+At).  Note  that  the  array  dimensions  used  in  amodel.F  are 
more  general  (m,n,kk)  since  they  can  represent  the  dimensions  of  any  mesh. 

Before  the  prognostic  tendencies  are  computed,  several  preliminary  diagnostic 
routines  are  called  to  define  arrays  used  later.  These  routines  include:  as/gdf.F which 
computes  sigma-coordinate  vertical  velocities  (sigdt),  ac//V3d.F  which  computes  the  two 
or  three-dimensional  divergence  (div3),  cgfemi.F  which  computes  and  adds  in  counter 
gradient  flux  terms  to  the  potential  temperature  and  water  vapor  mixing  ratio  when 
iamxgi=4  or  5,  and  ad/agf/r.F  which  computes  Smagorinsky  type  horizontal  eddy  mixing 
coefficients  when  iahsgm=1  (eks).  The  surface  parametertization  routine,  sfcflx.F,  is 
also  called  in  advance  to  produce  values  for  the  surface  fluxes  (hfluxi,  hfluxs)  and 
stresses  (stresx,  stresy)  used  as  lower  boundary  conditions  for  the  water  vapor, 
potential  temperature  and  wind  components  respectively. 

Next,  the  predicted  future  values  of  the  time  dependent  variables  are  obtained. 
The  terms  associated  with  the  slow  or  advective  modes  in  COAMPS  are  computed  first 
in  subroutines  afore.F{e3),  aforth.F (thrhs),  aforqx. F (qcZ,  qi3,  qr3,  qs3,  qv3),  arhsu.F 
(urhs),  arhsv. F(vrhs),  arhsw. F (wrhs),  and  ar/7sp. F(prhs).  In  the  above  routines, 
time  level  3  arrays  for  TKE  and  moisture  are  updated  because  they  do  not  contain 
terms  related  to  the  faster  moving  sound  and  gravity  waves  which  are  computed 
separately.  The  forcing  associated  with  the  fast  modes  are  integrated  in  subroutine 
a//7s.F  using  a  time-splitting  and  semi-implicit  computation  over  a  shorter  time  step 
determined  from  namelist  variable  mtaua.  Here  the  time  level  3  arrays  for  the  potential 
temperature,  wind  components  and  pressure  are  updated  adding  in  the  previously 
computed  advective  mode  tendencies  (thrhs,  urhs,  vrhs,  wrhs  and  prhs). 

Once  the  time  level  3  arrays  have  been  updated,  adjustments  are  made  to  the 
prognostic  variables.  These  adjustments  include  explicit  moist  physics,  radiation, 
turbulent  vertical  mixing,  and  temporal  filtering.  For  example,  subroutine  kfpara.F 
comprises  a  convection  parameterization  scheme  used  when  the  logical  variable 
lcupar=.true.  and  when  the  horizontal  resolution  of  the  nest  is  greater  than  that 
specified  by  namelist  parameter  dxmeso.  Subroutine  adjtq.F  ls  a  cloud  microphysical 
parameterization  scheme  that  handles  subgrid  scale  moisture  processes  when 
imoist=,true..  Subroutine  rad/af. F  parameterizes  the  long  and  short  wave  radiation 
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effects  upon  the  potential  temperature  when  logical  variable  Irad=.true..  Also  listed  as 
a  parameterization  subroutine  is  the  TKE  prediction  routine  afore.F  called  when 
ltke=.true..  This  routine  consists  of  1 .5  order  boundary  layer  closure  With  several 
options  for  parameterizing  the  eddy  mixing  coefficients  and  the  turbulent  mixing  length. 

Further  adjustments  are  made  to  the  time  level  3  arrays  by  applying  boundary 
conditions,  in  subroutines  abdye.F,  abdytq.F,  abdyt.F,  abdyuv.F,  abdyw.F,  abdyp.F, 
and  performing  vertical  mixing,  in  subroutines  amixtq.F,  amixnf.F,  amixuv.F,  amixw.F 
(ltke=.true.).  Finally,  temporal  oscillations  associated  with  the  leap  frog  integration 
scheme  are  smoothed  in  subroutine  atimef.F.  To  maintain  numerical  stability,  the 
damping,  diffusion,  filtering  and  horizontal  mixing  computations  are  performed  on  the 
time  level  1  arrays,  and  vertical  mixing  is  done  implicitly  on  the  time  level  3  arrays.  All 
other  quantities  are  computed  on  time  level  2  arrays.  See  Section  2  for  a  list  of 
references  that  describe  the  parameterization  and  adjustment  schemes  used  in 
COAMPS,  and  Appendix  D  for  further  descriptions  of  the  namelist  input  parameters. 

The  coding  in  the  above  subroutines  completes  one  model  iteration  with  updated 
predicted  values  stored  in  the  time  level  3  arrays.  Control  returns  to  the  forecast  routine 
where  the  atmospheric  model  is  called  for  up  to  six  additional  inner  nests,  thus 
repeating  the  amode/.F coding  sequence.  At  the  end  of  each  coarse  mesh  (1)  iteration, 
coamm.F  queries  whether  the  forecast  fields  are  written  to  ‘restart  files  (in  subroutine 
iomdI.F)  and/or  graphical  display  files  (in  subroutine  output.F)  as  described  in  the  next 
section. 

3.3.2.3.  Writing  Output 

COAMPS  produces  several  forms  of  model  output  in  subroutine  output.F  for 
viewing  case  study  results  and  testing  code  development.  Output  options  include:  2D 
horizontal  and  vertical  slices  of  data  viewed  in  numerical  format,  called  ‘quick  prints’,  as 
well  as  1 D,  2D,  and  3D  data  post-processed  through  a  variety  of  separate  graphics 
programs.  Figure  1 1  depicts  the  flow  diagram  for  writing  out  COAMPS  results.  Here, 
the  subroutines  are  briefly  introduced  beginning  with  an  overview  of  2D  data  output. 
More  detailed  information  concerning  COAMPS  graphics  is  presented  in  a  separate 
document. 

In  subroutine  aprint.F,  a  predetermined  set  of  model  fields  is  written  to  standard 
output  in  2D  numerical  quick  print  format  at  time  intervals  specified  by  the  user 
(kprnta).  The  standard  output  file  also  contains  standard  error  messages,  and  thus,  in 
addition  to  numerical  output  of  model  fields,  also  indicates  if  the  job  has  completed 
successfully.  The  2D  horizontal  and  vertical  slice  data,  generated  for  graphical  display 
purposes,  are  specified  by  the  user  before  a  COAMPS  model  run.  This  information  is 
read  into  COAMPS  through  graphics  ‘directives’  given  in  files  whose  names  correspond 
to  the  input  parameters  npfil  (horizontal  slice  data)  and  xsfil  (vertical  slice  data). 
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outputf 

write  forecast  data 


enter 


N^prntaSY 


ocards 

iteration^ 


xcards^'^ 

“\teration/ 


"N^ratio^/ 


aprinth  write  2D  horizontal 
and  vertical  slice  “quick  print" 
data  to  standard  output 


outsfc.f,  aoutp.f,  aoutz.f:  write 
2D  horizontal  slice  fields;  surface, 
constant  height  or  pressure  levels  to 
files  for  plotting  in  pbtxy 


aoutxz.t  write  2D 
vertical  slice  fields  to  files 
for  plotting  in  coamx 


asavid.f:  write  1D  grid 
point  data  to  ‘tee/^  files  for 
plotting  in  coamid 


return 


Figure  1 1 :  Flow  chart  of  the  subroutines  called  by  coa/nm./ that  write  out  model  data.  The 
input  parameters  kpmta,  tld,  ksaves,  tvis  represent  forecast  times,  defined  in  Appendix 
D,  that  are  converted  into  iteration  numbers,.  The  ocards  and  xcards  iterations  refer  to 
times  given  in  the  2D  graphics  directive  files  created  by  the  user. 
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The  2D  graphics  directives  are  used  in  subroutines  outsfc.F,  aoutp.F,  aoutz.F,  and 
aoutxz.F,  to  produce  the  following  types  of  output; 

•  outsfc.F-  surface  fields 

•  aoutp.F  -  horizontal  slices  of  fields  at  pressure  levels 

•  aoutz.F  -  horizontal  slices  of  fields  at  height  levels 

•  aoutxz.F-  vertical  slices  of  fields 


The  output  files  generated  are  named  following  a  standardized  naming  convention. 
Each  unique  filename  uses  36  characters  to  identify  1 1  pieces  of  information  about  the 
data  it  contains.  For  example,  filename  wspda21 99406081 200600000002000000hsJ 
represents: 


1994060812 

0060000 

00020 

00000  hs  1, 

L  - 3 

■ 

type  of 
level 

date/time 

group 

forecast 

time 

l“and2“ 
height  level 

type  of 
slice 


The  filename  indicates  that  the  data  are  of  horizontal  wind  speed  from  the  atmospheric 
model  on  the  second  nest.  Furthermore,  the  data  are  a  horizontal  slice  on  a  height 
surface  at  a  height  of  20  m,  6  hours  into  a  forecast  initialized  at  12Z  on  June  8th,  1994. 
Table  4  summarizes  the  information  given  by  the  output  filename. 


Table  4 :  COAMPS  Output  Fi' 

.e  Naming  Convention 

.OUTPUT  FILENAME 

TYPE 

POSITION 

OPTIONS-DESCRIPTION 

field  name  indicator 

a4 

1-4 

see  COAMPS  graphics  document 

model  type  indicator 

a1 

5 

a  -  atmosphere 

0  -  ocean 

nest  number  indicator 

il.1 

6 

1  -  coarse  mesh 

2,3, 4, 5,6  or  7  -  inner  nests 

alO 

7-16 

YYYYMMDDHH 

forecast  time:  hr 

i3.3 

17-19 

a  time  within  the  forecast  period 

min 

12.2 

20-21 

evenly  divisible  by  the  time  step 

sec 

i2.2 

22-23 

level:  meters 

i5.5 

24-28 

height  level  (horizontal  slice  data) 

second  level:  meters 

i5.5 

29-33 

2"°  height  level  (for  thickness  fields) 

level  type  indicator 

a2 

34-35 

si  -  mean  sea  level 
hs  -  height  surface 
ml  -  marine  level 
sf  -  surface 
pr  -  pressure  level 

plot  type  indicator 

a1 

36 

1  -  horizontal  slice 

X  -  vertical  slice 
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Subroutines  for  producing  1 D  (grid  point)  and  3D  volumes  of  model  data,  are 
briefly  mentioned  here.  The  1D  data  output  is  determined  from  the  coamni  namelist 
variables  tid,  dsetid,  cpar,  city,  ild  and  j1d.  Subroutine  asavId.F creates  1D 
graphics  files  (with  filename  prefix  ‘tsei’)  at  time  intervals  given  by  tid.  And  subroutine 
iosig.F creates  3D  volumes  of  sigma  level  data  (with  the  filename  ending  in  ‘sgf)  and 
2D  surface  data  (with  the  filename  ending  in  ‘sff)  at  time  intervals  given  by  two  namelist 
parameters:  icycle  (for  use  in  data  assimilation)  and  ksaves  (for  graphical  display: 
Viz5d,  GrADS).  The  standard  set  of  output  produced  by  this  subroutine  is  listed  in 
Table  5  for  reference. 

Table  5 :  Sigma  level  data  and  surface  fields  written  out  by 


COAMPS 


Volumes  r  .  .  .... 

2D  Surface  Fields 

pressure  (Pascals) 
perturbation  pressure  (Pascals) 
potential  temperature  (K) 
u-momentum  component  (m/s) 
v-momentum  component  (m/s) 
water  vapor  mixing  ratio  (kg/kg) 
cloud  water  mixing  ratio  (kg/kg) 
ice  crystal  mixing  ratio  (kg/kg) 
rain  drop  mixing  ratio  (kg/kg) 
snow  flake  mixing  ratio  (kg/kg) 
turbulent  kinetic  energy  (m**2/s**2) 
vertical  velocity  (m/s) 

ground  wetness  (fraction) 
ground/sea  surface  temperature  (K) 
boundary  layer  height  (m) 
sensible  heat  flux  (W/m**2) 
latent  heat  flux  (W/m**2) 
wind  stress  (Newton/m**2) 

Another  option  is  available  for  viewing  primarily  2D  model  results  using  the 
graphics  package  Visual.  Subroutine  visout.F creates  VisualtWes  (with  the  filename 
ending  in  ‘GRDf,  GRD2,  or  GRD3)  at  time  intervals  given  by  namelist  inputs  tvisi, 
tvis2. 


This  completes  the  overview  of  the  COAMPS  code.  At  this  point,  the  flow 
structure  of  the  code,  the  subroutine  functions,  and  the  use  of  certain  namelist 
parameters  have  been  described.  In  the  next  section,  the  procedure  for  submitting  a 
COAMPS  job  is  reviewed. 
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4.  Execution  Procedure 


The  execution  procedure  for  a  COAMPS  job  involves  two  phases:  case  study 
setup  and  submitting  a  COAMPS  job.  Before  a  COAMPS  job  begins,  several 
preliminary  steps  are  taken  involving  preparation  of  the  initial  conditions  and  the  code. 
The  steps  outlined  below  show  how  to  get  started  using  COAMPS  as  a  new  user,  and 
how  to  set  up  a  new  case  study.  First  time  users  of  COAMPS,  need  to  obtain  the 
COAMPS  tar  files,  set  up  the  COAMPS  directory  structure  and  compile  the  master 
COAMPS  libraries  and  executables.  These  steps  are  covered  in  Section  3.2.  Also 
described  in  Section  3.2  are  the  steps  for  obtaining  two  templates:  a  run  script  used  to 
run  the  COAMPS  executables,  run.coamps,  and  a  makefile  used  to  create  modified 
executables,  Makefile.  The  steps  given  here  describe  how  to  run  a  user-customized 
COAMPS  simulation  with  the  option  of  modifying  the  original  source  code. 

4.1 .  Preparation  Steps 

CODE  MODIFICATIONS 

(To  run  unmodified  source  code,  goto  step  10.) 

1 .  cd  to  the  subdirectory  /coamps/mod 

2.  type  “mkdir  modsd’  where  modsd  ls  any  user-defined  subdirectory  name  that  will 
contain  user  generated  code  modifications 

3.  cd  to  this  new  subdirectory  modsd 

4.  type  “cp  .yremplates/Makefile to  get  the  Makefile  template  into  the  ‘modsd 
subdirectory. 

5.  copy  subroutines  from  any  of  the  subdirectories  in  /coamps#/Hbsrc  for  modification 

6.  modify  and  add  subroutines 

7.  copy  all  include  files  {*.h)  from  the  subdirectories  in  /coamps#/libsrc that  are  called 
by  the  subroutines  that  have  been  modified 

8.  edit  the  Makefile  to  include  object  files  corresponding  to  subroutines  that  have  been 
modified  or  added 

9.  make  the  new  executable:  type  “make  machine  name”  (See  Table  3) 

INITIAL  INPUT  FIELDS 

1 0.  obtain  the  initial  input  fields 

•  <real  data  case> 

To  perform  a  NOGAPS  ‘cold  start’,  COAMPS  uses  1°  resolution  NOGAPS 
fields  In  flat  file  format  (32-bit,  IEEE)  and  ADP  observational  data  in  FGGE 
format.  These  data  are  available  from  FNMOC  and  can  be  obtained  using  an 
anonymous  ftp.  For  a  typical  model  run,  a  NOGAPS  analysis  or  forecast 
provides  the  first  guess  fields  and  12  hourly  NOGAPS  forecasts  provide 
boundary  conditions  throughout  the  duration  of  the  simulation  study  period. 

In  the  absence  of  NOGAPS  forecasts,  fixed  boundary  conditions  are  used 
(although  this  is  not  recommended).  Once  these  data  are  obtained,  they 
should  be  copied  into  the  appropriate  subdirectory: 
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/coampsWdatabase/nogaps  and  /coamps#/database/adp.  Other  databases 
used  in  COAMPS  (surface  fields,  terrain  characteristics  and  climatology  data) 
are  provided  with  the  initial  COAMPS  tar  file  database.tar.  These  data  are 
also  required  to  be  in  flat  file  format.  (Note  that  the  subdirectory  locations  of 
this  initial  input  data  must  be  provided  as  namelist  input:  parameters  dsngff, 
dsclim,  dsdted,  dsgiss.) 

•  <idealized  case> 

Obtain  or  create  an  idealized  sounding  (pressure,  temperature,  moisture  and 
winds)  and  input  these  data  through  namelist  arrays  psnd,  tsnd,  qsnd  usnd, 
vsnd.  Perturbations,  used  to  modify  these  initial  idealized  profiles,  may  be 
defined  by  the  user  in  subroutine  flds##.F .  Perturbations  are  added  when 
the  namelist  parameter  icase=##  coincides  with  the  number  in  the  flds##.F 
filename.  The  surface  fields  can  be  defined  from  either  user-specified 
namelist  parameters:  aindpct,  seatmp,  zrough,  albdo,  sfcwet  or  from  the 
standard  COAMPS  surface  databases  located  in  subdirectory 
/coampsWdatabase.  (Note  that  if  real  surface  data  is  used,  the  subdirectory 
locations  of  the  surface  databases  must  be  provided  as  namelist  input: 
parameters  dsclim,  dsdted,  dsgiss.) 

CASE  STUDY 

11  .cd  to  the  subdirectory  /coamps/run 

12.  type  “‘mkdir  runsd’  where  runsd  'is  any  user-defined  directory  name  that  will  contain 
a  case  study  run  script  and  output  files 

13. cd  to  this  new  subdirectory  runsd 

14.  type  ‘cp  .yremplates/run.coamps  to  get  the  COAMPS  run  script  template  into 
the  mnsc/ subdirectory. 

15.  edit  the  run.coamps  run  script  to  customize  the  namelist  input  for  a  particular  case 
study  (i.e.  modify  the  gridnl,  coamni,  and  mvoini  namelists  (See  Section  4.2)) 

16.  if  the  COAMPS  executables  have  been  modified,  make  sure  that  the  binDir 
variable  in  the  run.coamps  script  points  to  the  location  of  the  new  executables. 

17.  create  a  /coamps/data  subdirectory  for  graphical  and  model  output  data 

These  steps  summarize  the  preparation  phase  of  a  COAMPS  job.  Correct 
specification  of  user  input  (step  1 5)  is  of  primary  importance  in  the  above  sequence. 
Because  COAMPS  includes  a  variety  of  important  features  and  options  as  namelist 
input,  the  user  has  the  ability  to  interact  with  the  code  through  these  namelists. 

Although  this  is  of  great  benefit,  high  levels  of  user  input  also  creates  the  opportunity  for 
error.  To  educate  the  user  and  reduce  the  risk  of  error,  categories  of  user  input  are 
grouped  below  based  upon  their  function  with  the  COAMPS  system. 
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4.2.  COAMPS  Input  File:  Gridnl,  Coamn!  and  Mi^o/n/ Namelists 

Perhaps  the  single  most  important  step  in  running  COAMPS  successfully  is 
indicating  the  proper  namelist  input.  Many  of  the  more  important  model  parameters  are 
specified  by  the  user  in  the  three  model  namelists  gridnl,  coamn!  and  mvoinl.  The 
parameters  from  these  namelist  and  their  user-specified  values  are  assigned  in  a  single 
input  file  which,  for  purposes  of  discussion,  is  called  nlist.input.  This  input  file  is  created 
when  the  COAMPS  run  script  is  submitted  (See  Appendix  B  for  a  sample  input  file).  A 
complete  list  of  the  namelist  input  parameters,  default  values,  and  definitions  is  given  in 
Appendices  C  and  D  arranged  alphabetically  according  to  the  type  of  input  (i.e.  integer, 
real,  character,  logical).  Categories  of  user  inputs  described  here  include  the  following 
parameter  groups:  model  setup,  case  study  setup,  and  data  manipulation. 

4.2.1.  Model  Setup 

This  section  describes  the  inputs  that  allow  the  user  to  set  up  the  model  domain 
and  grid  structure,  the  integration  procedure,  and  the  boundary  conditions  for  a 
particular  case  study. 

4.2.1 .1.  Grid 

Specification  of  the  namelist  parameters  associated  with  the  grid  setup  are 
critical  In  producing  a  proper  model  run.  Their  values  vary  from  case  to  case  depending 
on  the  region/domain  of  interest,  and  the  vertica:l/horizontal  spacing  required  to  resolve 
meteorological  features  of  interest. 

Input  parameters: _ 

the  number  of  nests  nnest;  horizontal  grid  spacing  of  the  coarse  mesh  deix,  dely; 
vertical  grid  spacing  dsigma;  number  of  grid  points  for  each  nest  ma(nn),  na(nn),  kka; 
map  projection  nproj;  standard  latitude(s)  and  longitude  phnti ,  phnt2,  ainnt;  reference 
latitude  and  longitude  riat,  riong;  reference  grid  points  iref(1),  jref(1);  reference  grid 
points  of  inner  nests  ii(nn),  jj(nn)  (where  nn  is  the  nest  number) _ 


It  is  possible  to  determining  the  model  domain  in  advance  of  a  COAMPS  job  by 
running  the  script  file  domain.sh  located  in  the  /utility/T emplates  subdirectory.  Edit  this 
file,  inputting  values  for  the  grid  setup  parameters,  and  then  type:  “domain.sh”.  A  meta 
file  of  the  nested  model  domain  is  created  and  displayed  on  the  screen,  provided  that 
the  DISPLAY  environmental  variable  has  been  specified  and  the  2D  graphics 
executable  plotxy  is  available.  Once  you  are  satisfied  with  the  model  domain,  insert 
those  parameter  values  into  the  coamn!  namelist.  Examples  of  three  model  domains 
are  given  in  Appendix  E.  Tables  E-1  ,E-2,E-3  give  the  namelist  parameter  values  used 
to  set  up  a  model  grid  based  upon  the  mercator,  lambert  conformal,  and  polar 
stereographic  map  projections  respectively.  Figures  E-2,E-4,E-6  show  the 
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corresponding  COAMPS  model  domains  obtained  by  running  domain.sh  script  for  each 
of  the  three  examples. 

4.2.1. 2.  Integration 

COAMPS  uses  second-order  leapfrog  time  differencing  to  integrate  the  slow 
modes  (advective  solution)  using  a  time  step,  delta,  that  depends  upon  the  grid 
resolution  and  advective  speeds  to  ensure  computational  stability.  Because  COAMPS 
is  nonhydrostatic,  the  simulations  also  contain  sound  and  gravity  wave  solutions.  Their 
faster  propagation  rates  require  integration  over  a  smaller  forward  time  step,  designated 
by  delta/mtaua  where  mtaua  is  the  number  of  small  time  steps  per  large  time  step. 

Two  important  numerical  modeling  improvements  allow  for  larger  values  of  the  small 
time  step.  These  include  semi-implicit  vertical  differencing  and  the  quasi-compressible 
assumption  which  permits  the  speed  of  sound  to  be  much  less  than  its  typical 
atmospheric  value  of  -350  m/s. 

Input  parameters: _ 

the  starting  and  ending  forecast  times  ktaust,  ictauf;  the  large  time  step  delta;  the 
number  of  small  time  steps  per  large  time  steps  mtaua;  and  the  speed  of  sound  cbar 


4.2.1. 3.  Boundary  Conditions 

These  parameters  allow  for  specification  of  the  lateral  boundaries  and  damping 
at  the  upper  boundary.  Five  types  of  lateral  boundary  conditions  are  available  for 
idealized  simulations  and  two  types  are  available  for  real  data  simulations  each  defined 
in  Appendix  D.  Real  data  boundaries  can  be  computed  by  the  Perkey-Krietzberg 
method  (ibdya=jbdya=nbdya=6),  which  specifies  the  boundaries  according  to  the 
NOGAPS  tendencies,  or  by  the  Davies  method  (ibdya=jbdya=nbdya=7),  which  uses 
the  actual  NOGAPS  fields.  In  general,  the  Davies  method  is  recommended  for  most 
real  data  model  runs.  Additionally,  two  types  of  upper  level  sponge  layers  have  also 
been  developed  to  reduce  the  amplification  of  vertical  wave  reflections  off  the  model 
top. 

Input  parameters: _ 

damping  upper  layer  Iralee,  Ispong,  nrdamp,  rdtime;  lateral  boundary  condition 
indicators  ibdya,  jbdya,  nbdya;  gravity  wave  speed  (for  idealized  radiative  horizontal 
boundary  conditions)  cstar;  number  of  grid  points  in  real  data  lateral  boundary 
computations  nbdypt;  times  for  obtaining  coarse  mesh  (1)  lateral  boundary  tendencies 
from  NOGAPS  fields  itauin,  itauf;  time  interval  for  updating  coarse  mesh  (1)  lateral 
boundaries  during  COAMPS  forecast  kgetbc _ 
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4.2.2.  Case  Study  Setup 


These  parameters  govern  the  model  features  that  are  active  for  a  particular  case 
study.  They  are  grouped  into  three  categories:  analysis  fields  parameters,  model 
physics  parameters  and  numerical  damping  parameters.  This  group  of  inputs 
determine  the  following  characteristics  of  the  case  study:  type  of  model  run  (idealized  or 
real  data  run),  data  assimilation  features  (type  of  update,  MVOI  analysis),  physical 
parameterizations  active  (TKE,  radiation,  moisture,  surface  fluxes,  surface  energy 
budget),  numerical  damping  options  (diffusion,  time  filter,  divergence  damping,  upper 
level  sponge  layers),  in  addition  to  miscellaneous  dynamical  features  (advection, 
geostrophic  wind).  Note  that  these  inputs  can  lead  to  complications  when  poorly 
specified. 

4.2.2.1.  Analysis  Fields 

These  parameters  allow  for  the  specification  of  fields  required  to  run  the  model. 
Idealized  initial  conditions  are  more  simply  specified  than  are  real  data  initial  conditions 
because  they  are  set  up  by  the  user  in  the  COAMPS  input  file.  Perturbations  from  a 
user  generated  reference  sounding  are  fed  into  the  COAMPS  analysis  routine 
(coama.F)  through  a  f/ds##.Ffile.  In  contrast,  real  data  setup  is  considerably  more 
complex  using  a  variety  of  sources  to  prepare  the  initial  conditions  including  NOGAPS 
analysis  or  forecasts,  COAMPS  history  fields,  climatology  fields,  surface  parameters, 
and  observational  data. 

Input  parameters: _ _ 

For  idealized  model  runs:  the  case  study  indicator  icase^^O;  latitude  of  f-plane  flat; 
interpolation  flag  for  nested  fields  jcm2fg;  idealized  surface  conditions  iindfig,  iseaflg, 
izOfIg,  iaibfig,  igwtfig,  aindpct,  seatmp,  zrough,  albdo,  sfcwet,  itopofig;  idealized 
reference  sounding  umean,  vmean,  ipsnd,  itsnd,  iqsnd,  iuvsnd,  psnd(k),  tsnd(k), 
qsnd(k),  usnd(k),  vsnd(k)  (where  k  is  the  vertical  level  index); 

For  real  data  model  runs:  the  case  study  indicator  lcase=0;  type  of  update  iupd;  multi¬ 
variate  optimum  interpolation  analysis  (MVOI)  loi,  lolmf;  nonlinear  variational  balancing 
ivar;  directory  locations  of  initial  input  dsngff,  dsnrff,  datfil,  dsclim,  dsdted,  dsgiss; 
real  data  surface  conditions  Inrhty,  Itopoa,  nftopo,  gwetmx,  zOzOmx;  data 
assimilation  icycle,  itauin,  itauf,  kgetbc,  loi,  iupd:  Cressman  analyses  ItanI,  IqanI, 
Itanis 


Note  that  real  data  surface  conditions  may  also  be  used  for  idealized  model  runs,  and 
idealized  surface  conditions  may  be  used  for  real  data  runs. 
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4.2.2.2.  Model  Physics 


COAMPS  contains  many  inputs  that  govern  various  parameterization  options  for 
specifying  the  surface  and  boundary  layers,  moist  physics,  convection  and  radiation 
schemes.  Refer  to  Appendix  D  for  further  definitions, 

input  parameters: _ _ 

surface  layer  processes  Iflux,  Isfcen,  lies;  boundary  layer  processes  Itke,  iahsgm, 
iamxgl,  iashsm,  shO,  smO,  lies;  moist  and  convective  processes  iisp,  Imoist,  lice, 
Icupar,  cupd;  radiative  processes  Irad,  dtrad,  njump,  Icool;  vertical  advection  of  TKE, 
rain  and  snow  dxmeso,  dzqlim _ 


4.2.2.3.  Numerical  Damping 

To  control  numerical  instability  and  the  growth  of  unrealistic  tendencies,  several 
damping  features  have  been  added  to  the  COAMPS  code. 

Input  parameters: _ 

diffusion  and  filtering  Idiff,  dif2nd,  dif4th,  Ifft;  divergence  damping  Iddamp,  divdmp; 
time  filtering  robert;  upper  level  sponge  layer  Iralee,  rdtime,  nrdamp,  Ispong _ 


4.2.2.4.  Additional  Model  Features 

These  parameters  provide  options  for  the  advective  scheme  (e.g.  2"'^  order,  4**^ 
order,  flux  or  advective  schemes)  and  allow  for  adding  a  large  scale  pressure  gradient 
via  the  geostrophic  wind. 

Input  parameters: _ 

advection  scheme  iadvct;  geostrophic  wind  Ivgeo,  ugeo,  vgeo _ 


4.2.3.  Data  Manipulation 

Several  namelist  parameters  are  used  to  control  the  flow  of  data  into  and  out  of 
the  COAMPS  system.  These  inputs  are  categorized  into  three  groups:  data  input  and 
output  parameters  (used  for  analysis,  restarts  and  data  assimilation  purposes),  graphics 
parameters  (used  to  view  model  output),  and  diagnostic  parameters  (used  to 
numerically  evaluate  model  results). 

Input  parameters: _ 

directory  locations  dsngff,  dsnrff;  times  for  printing/saving  output  kprnta,  ksavea, 
ksaves,  isavefrq,  ifsave;  data  format  type  Iwritu;  1 D  graphics  output  dsetid,  tid, 
i1d(nn),  j1d(nn),  cpar(nn),  city(nn);  2D  graphics  output  npfil,  xsfll,  Iwritu;  3D  graphics 
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output  tvisi,  tvis2,  vpar1,  vpar2,  Iprvis,  prsvis,  hgtvis,  ksaves,  Iwritu;  diagnostic 
output  laper,  lasym,  Idgadj,  Icuppr,  [prnti _ 


Although  most  of  these  output  times  are  designated  by  time  intervals,  note  that 
parameter  ksaves  is  an  array  containing  actual  times  during  the  forecast  that  3D 
volumes  of  sigma  level  data  are  to  be  saved.  Caution  is  advised  when  outputting  3D 
volumes  over  a  long  forecast  period  or  at  many  times,  especially  for  more  than  one 
mesh  since  they  are  quite  large  (a  (65x65x30)  grid  requires  approximately  100Mb  of 
space). 

Having  reviewed  the  more  important  COAMPS  input  parameters  in  this  section, 
and  using  Appendix  B  as  an  example  of  a  COAMPS  input  file,  the  user  is  now  able  to 
create  a  COAMPS  input  file  that  sets  up  any  type  of  model  run.  The  final  step  in 
performing  a  COAMPS  simulation  is  to  submit  the  run  job  as  described  in  the  next 
section. 

4.3.  Submitting  a  COAMPS  Job 

The  user  is  ready  to  submit  a  COAMPS  job  once  the  following  steps  have  been 
completed: 

1 .  The  analysis  and  forecast  codes  have  been  compiled  and  linked  to  form  the 
COAMPS  executables  {coamps_analysis.exe  and  coamps_forecast.exe. 

2.  The  data  has  been  obtained  to  specify  the  Initial  input  fields. 

3.  The  COAMPS  run  script  run.coamps  has  been  edited  to  update  the  namelist  input 
parameters  (i.e.  the  namelist  input  file  nlistinput  has  been  modified  to  suit  the 
particular  process  study). 

In  the  example  given  below,  the  COAMPS  run  script  template,  run.coamps,  has  been 
modified  and  renamed  mn.mnsc/ where  runsdls  a  user-created  subdirectory  under /mn. 
This  script  creates  a  COAMPS  namelist  input  file  called  nlist.$ddtg.runsd  where  $ddtg 
is  the  a  10  digit  date-time  group  (YYYYMMDDHH)  provided  to  the  script  file  as  an  input 
argument. 

Within  the  /run/mnsd  subdirectory,  submit  the  run  script  by  typing:  “run.runsd 
$ddtg  3”.  The  first  input  argument  represents  the  starting  date-time  group.  The 
second  is  a  run  action  argument  indicating  whether  the  analysis  or  forecast  is  to  be  run. 

Table  6 :  Options  for  the  2nd  run  script  input  argiiment 


ARGUMENT 

ACTION 

0 

no  action 

1 

run  analysis 

2 

run  forecast 

3 

run  analysis  and  forecast 
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When  the  second  input  argument  is  ‘1’  or  ‘3’,  the  run  script  executes  the  analysis 
program  using  the  following  command: 

$HOME/coamps/bin/coamps_analysis.exe  <  nlist.$ddtg.runsd  >  output.a.$ddtg 

The  user  input  is  redirected  into  the  analysis  program  from  file  nlist$ddtg.runsd  and 
standard  output/error  is  redirected  out  to  file  output.a.$ddtg.  Upon  completion,  the 
directory  specified  by  namelist  input  dsnrff  contains  the  analysis  fields  produced  by  the 
executable.  Scan  this  directory  and  view  the  file  output.a.$ddtg  to  ensure  that  the 
analysis  completed  successfully. 

Similarly,  when  the  second  input  argument  is  ‘2’  or  ‘3’,  the  run  script  executes  the 
forecast  using  the  following  command: 

$HOME/coamps/bin/coamps_forecast.exe  <  niist.$ddtg.runsd>  output.m.$ddtg 

User  input  is  redirected  into  the  forecast  program  from  file  nlist.$ddtg.runsd  and 
standard  output/error  is  redirected  out  to  file  output.m.$ddtg.  Upon  successful 
completion  of  the  COAMPS  forecast,  user-requested  forecast  fields  reside  in  the 
directories  specified  as  namelist  input  dsnrff  and  dsetid.  These  data  may  be  viewed 
graphically  with  a  variety  of  display  packages  that  process  1 D,  2D  or  3D  model  output 
the  details  of  which  are  covered  in  a  separate  document.  For  a  quick  analysis  of  the 
COAMPS  results  at  time  intervals  given  by  the  user  input  parameter  kprnta,  the  user 
may  look  at  ‘quick  prints’  of  horizontal  and  vertical  slice  numerical  model  output 
embedded  in  the  standard  output/error  file  outputm.$ddtg.  The  output  to  this  file  also 
indicates  whether  the  job  has  completed  successfully. 

A  data  assimilation  update  cycle  may  be  submitted  by  inputting  two  additional 
run  script  arguments.  The  third  input  argument  is  a  1 0  digit  ending  date-time  group 
(YYYYMMDDHH)  and  the  forth  argument  is  the  number  of  hours  between  updates 
(typically  12  hours).  To  submit  a  12  hour  data  assimilation  cycle,  the  user  types: 
“run.runsd  $ddtg  3  $enddtg  12”.  The  command  initiates  a  forecast  starting  at  $ddtg 
and  continuing  through  $enddtg  assimilating  data  every  12  hours.  This  type  of  run  is 
only  appropriate  when  observations  are  available  every  12  hours  throughout  the  time 
period  since  these  data  are  needed  to  adjust  the  first  guess  fields. 
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Appendix  A:  Sample  makefile 


This  Makefile  template  is  from  the  /coamps#/mod/remplates  directory  that  is  created  by 
running  the  script  file  get.templates  for  the  SGI  machine  Rossby. 

nic  MukeRh 


tV  makci'iic  lo  biiiid  pioLMiiiii  i;()aniT)s_anal_v.si.s.cxc  and  a<)anips_fori.'ca.st.c\i: 


Oilier  dependeneies 


iiieliidc  /diua..ro.ssb_\  6/usoiS('coaitips/coaiiips2.0/c()nliii.os 

include /data_rossby6/users/coamps/coamps2.0/config.user  ,  ^  ' 

LIBDIR  =  /data_rossby6/iisers/coamps/lib/2.0 

MAINANL  =/dala_rossby6/uscrs/coamps/coaraps2.0/src/coamps_analysis/coamps_analysis.F 
MAINFCST  =  /daia_rossby6/uscrs/coamps/coamps2.0/src/coamps_forecast/coamps_forecast.F 
INCL_PROL_MAIN  =/data_rossby6/users/coamps/coamps2.0/prologues/main 
INCL_PRQLjCOAMPS  =  /data_rossby6/uscrs/coamps/coamps2.0/prologues/coampslib 
INCL_PROL^OI  =  /data_rossby6/users/coamps/coanips2.0/prologues/oilib 


LIBS  =S(Lini:)IR),'coampslib.a  \ 

$(LIllDIR:i/rn()clib.a  \  'V,; 

.S(LIllDIR)/oilib.a  \  ^ 

S(LIBD.IR)/eoainpi;lib.a  \ 

S(]  JBDlR)/nI.J)oqlib.a  \ 

S(lJnDIR)/fisiipaklib.a 

MAEEFILE  =  Makefile 

#  OBJS;  include  the  coamps  .o  files  that  the  user  modified 


PROGRAMS  =  coamps_analysis.exe  \  , 
coamps_forecast.exe 

target :  $_(PROGRAMS):$COBJS)  $CLIBS) 


(continned) 


A-1 


#  loader  diriKtive 


coamps_ana1ysis.exe  :  $(M4INANL)  $(OBJS) 
$(RM)$@ 

S(LD)  $(LDELAGS)  -o  $@  -I$(INCLJE?ROO 
$CEXTRALIBS)  ^ 

coamps_forecasLexe ;  $(MAI]SDFCST)  S(OBJS)  $(LIBS) 

SOUPBLAGS)  -o  $@  -ilBfCXJHEtOtJS 
S(r.XTRALIBS)  .  y 

:ray  $@  .A./bto7$;@„  .  : 


i$CFC)  -c  $(FFLAGSM)  -I$(INCL_PROL_COAMPS)  -I$(INqLJPROL_OI)  $*JF 


Tf  common  rules  for  ali  Makefiles  -  do  not  edit 

F rHfifT:Tr rrfriJ  "  rr  rff#  fi  4-rf^'f:  •rri^TIrr  iritith  t^?nnt  tf  iiiHrff  if  f^t fi  ff  ff  rr/l"l  ‘rmV’f?  fftrlf 

emptyrule: 


#  done  ■  /■' '  ■ 


AppsndiX  B:  sample  coamps  user  input  file  niist.input 


file  nlisLinput _ 

&gridnl 

kka  =  30, 
ma  =  (^,73,  109, 
na  =  65,  73,  91, 
iibdypt=  7, 
nnest  =  3 , 

&eiid 

&coanffil 


dsclim  = 'insert  directory  path  of  climatology  data', 

dsdted=  inscit  diiectoiy  jxuii  oi"  1km  terrain  data’, 

dsgiss  =  'insert  director^'  psith  of  Goddard  Institute  of  Space  Studies  data', 

dsngff  =  'insert  directory  path  of  NOGAPS  flat  files', 

dsnrfif  =  'insert  directory’  path  of  COAMPS  flat  files', 

npfil  =  'insert  name  of  ocards  file;  2D  horizontal  graphics  directives'  , 
xsfil  =  'insert  name  of  xcards  file;  2D  vertical  graphics  directives', 

icase  =  0, 

ibdya  -  7,  . 

jbdya  =  7*  ■ 

kgetbc=  12, 
iqjroj  =  2, 

alnnt  =  5,0, 

pMtl  =  60.0, 

phnt2  =  30.0, 

riatl  =  45.0, 

riongl=  5-0, 
iref  =  33, 
jfef  =  33?' 
ii  =21,20 
jj  18,17 

delx  =  SIOOO.O. 
dely=  81000.0. 
delta  =  240.0, 

ktaust=  0,0,0, 
ktauf  =  24, 0.  0. 

24,0.0. 

24,0,0, 

intaua  =  2,  .  -■■'■u. 

alpha  =  1.0, 

dsigma  =  7500.0,  5800.0,  4200.0,  2500.0,  1000.0, 

1000.0.  750.0,  750.0.  750.0.  750.0, 

:750.0,  750UO,  1000.0,  1000.0, 1000.0,  ? 

1000 A  800.0.  800.0,  800.0,  600.0, 

400:6,300.0,  200.0,140.0,  90.0, 

60.0,  40.0,  30.0,  20.0,  20.0, 
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(continued) 


Oj04, 

dief&=  -(XipOS, 
lddainp  =  t, 
divdmp  =  MtS, 
lralee  =  f, 

=  150.0, 

nrdamp  =  8, 

1X2, 

I2way  =  t 

lygep  =f, 
iaiivct=  2, 
cbar  =  240.0, 

itaoffl  =  12, 
itaiif  =  24, 
icycle=  12, 
tupd  =  0, 
loi  =:f, 
loimf  =f, 

Inrhfi'  =  f, 

Iqanl  =  f, 

Itahl-f, 
lpseud  =  f, 
ivar  =  0 

Icupar  =  t, 
lcupjar  =  f, 
Imoist=t, 
lice  =t, 

Irad  =t, 
njump=l, 

Icool  -f,  4 

Isfcen  =  t, 

Iflux  =t. 

Unit  =t, 

lies 

like  =t, 
almin  =  0.6^ 

alO  =  IJa 

shO  =  0.6;75, 

smO  —  04, 

iahsgm  =0, 
iamxgl=  4, 
iashsm=  4, 

ilndflg=l, 

ialbag  =  1, 

i20flg  =  1, 
igwtflg=  1, 
iseafig  =  1, 
itopoflg  =  1, 

Itopoa  =  t, 
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nftopo  =  3, 

kprnla  =  48, 0,  0, 
ksavea  =  144,  0,  0, 
ifsavc  =  2 

ksaves  =  24,  0,  0, 
hvritu  =  t, 

Ipmti  =  f, 
locard  =  f, 

Idbms  =f, 

Idbmsi  =  f, 

Idbmso  =  f, 

&end 

&rnvoinl 

datfil  =  'insert  directory’  path  of  observational  data’ 
&cnd 
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Appendix  C:  input  parameters:  grr/dh/ namelist 


VAR 

DEFAULT 

VALUE 

,  DEFINITION :  - 

INTEGER  VARIABLES 

kka 

30 

Number  of  vertical  height  levels  in  the  atmosphere 

NOTE:  (kka  must  be  equal  to  the  number  of  dsigma_ 
values  setup  in  the  namelist) 

Im 

16 

Number  of  vertical  pressure  levels  on  which  real  data 
initial  conditions  are  specified 

ma 

61 

ID  array:  Number  of  grid  points  in  the  x  direction  for 
each  mesh  NOTE:  (for  inner  nests  (ma-l)/3  must  be  a 
whole  number) 

na 

61 

ID  array:  Number  of  grid  points  in  the  y  direction  for 

each  mesh  NOTE:  (for  inner  nests  (ma-l)/3  must  be  a 
whole  number) 

nbdya 

7 

Real  data  lateral  boundary  condition  identifier  for  inner 
nests: 

(6)  Perkey-Krietzberg 

(7)  Davies 

nbdypt 

5 

Number  of  points  included  in  lateral  boundary 
computation  NOTE:  (ma  and  na  >  (2*nbdypt)) 

nnest 

1 

Total  number  of  meshes 
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AppGndiX  D:  input  parameters:  coamn/ namelist 


i  VAR 

DEFAULT  VALUE 

DEFINITION 

1  CHARACTER  VARIABLES 

city 

1  1 

Character*32, 2D  array:  City  names  associated  with 
grid  point  locations  (ild  jld)  for  each  mesh;  used  for 
labeling  ID  graphics  output  (getld.F,  asavld.F) 

cpar 

t  t 

Character*4, 2D  array:  ID  graphics  variable 
identification  string  for  each  mesh;  used  to  select 
variables  for  ID  graphics  output  {getld.F,  asavldF) 

dsclim 

'/home/hodur/masterr/' 

Character*80:  Directory  path  of  climatology  data; 
used  to  initialize  the  real  data  surface  parameters 
(albedo,  sst,  grdwet,  rland,  zsfc,  zO)  when  high 
resolution  data  not  available  (sfcpar.F,  soiltp.F) 

dsdted 

Yu/c/hodur/masdted/' 

Character*96:  Directory  path  of  1  km  terrain  data; 
used  for  horizontal  resolution  <  20  km  and  if  available 
(sfcpar.F) 

dsgiss 

Va/ops/etc/static/met 

/noraps/masgiss/' 

Character*80:  Directory  path  of  Goddard  Institute 

Space  Studies  high  resolution  albedo  and  zO  data 
(sfcpar.F) 

dsetng 

1  1 

Character*24:  ISIS  directory  path  of  NOGAPS  data 
input/output;  used  for  location  of  boundary  conditions 
and/or  cold  start  initial  conditions  when  Idbms=.true. 
(gethty.F,  getng.F,  getsst.F,  inhesc.F,  inngsc.F, 
innguv.F) 

dsetg 

1  t 

Character*24,  ID  array:  ISIS  directory  path  of 

CO  AMPS  data  input/output  for  each  mesh;  used  for 
location  of  COAMPS  files  when  ldbms=.trae. 

(gethty.F,  getng.F,  getsst.F,  inhesc.F,  inngsc.F, 
innguv.F) 

dsetld 

Character*80:  Directory  path  for  ID  graphics  output 
(getld.F,  asavldF) 

dsngff 

1  I 

Character*80:  Directory  path  of  NOGAPS  flat  files; 
used  for  location  of  boundary  conditions  and/or  cold 
start  initial  conditions  when  ldbms=.false.  (gethty.F, 
getng.F,  getsstF,  inngsc.F,  innguv.F) 

dsnrff 

1  1 

Character*80:  Directory  path  of  CO  AMPS  flat  files; 
used  for  location  of  CO  AMPS  input/output  needed  for 
incremental  update,  model  restart,  and  2D  and  3D 
graphics  output  when  ldbms=.false.  (coamm.F,  io*.F) 

gnamg 

1  1 

Character*32,  ID  array:  ISIS  geometry  name  for  each 
mesh;  used  by  ISIS  when  Idbms  and  ldbmsi=.trae. 

mditg 

t  1 

Character*32,  ID  array:  ISIS  model  type  for  each 
mesh;  used  by  ISIS  when  Idbms  and  ldbmsi=.tme. 

npfil 

'ocards' 

Character*80:  name  of  file  containing  2D  horizontal 
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VAR  DEFAULT  VALUE 


xsfil  'xcards' 


vparl 


vpar2 


REAL  VARIABLES 


albdo  0. 2 


alndpct  I  1.0 


climrh 


almin  I  0.02  m 


alnnt  240.0  ° 


100.0  m 


100.0  m/s 


.01,.01,.01,.02,.02,.10,.15 

,.20,.25,.30,.40,.50,.60,.70 

,.70,.70,84*.70 


300.0  m/s 


cupdx  10000.0  m 


delta  2.0  s 


DEFINITION 


level  plotting  information  (fields,  plot  type,  level,  time, 
mesh  number,  etc.) 


Character*80:  name  of  file  containing  2D  cross- 
sectional  plotting  information  (fields,  plot  type,  levels, 
time,  mesh  number,  etc.) 


Character^,  2D  array:  3D  graphics  variable 
identification  string  set  1  for  each  mesh;  used  to  select 
output  (visout.F) 


Character*4, 2D  arrays:  3D  graphics  variable 
identification  string  set  2  for  each  mesh;  used  to  select 
Visual  output  (visout.F) 


Idealized  value  of  surface  refractivity;  used  to  specify 
the  albedo  when  ialbflg  =  0  (user_rfce.F)  NOTE: 
(0.0  =  low  refractivity,  1.0  =  high  refractivity) 


Idealized  land/water  indicator;  used  to  specify  whether 
a  point  is  a  land  or  water  point  when  ilndflg  =  0 
(user_sfc.F)  NOTE:  (0.0  =  water,  1 .0  =  land) 


Minimum  mixing  length  factor  (afore.F) 


Longitude  aligned  north-south;  used  for  rotating  a 
lambert  conformal/polar  stereographic  projection 
(grid.F)  NOTES:  (increases  eastward  from  0°  to  360° 
beginning  at  Greenwich,  England)  (longitudes  outside 
the  model  domain  may  be  used  to  rotate  the  grid) 


Averaging  coefficient  for  semi-implicit  time 
differencing  (alhs.F,  alhsin.F). 


Mixing  length  constant/limiting  value  (coamm.F) 


Speed  of  sound  char=  .yjc^RT^  /  c„  (alhs.F,  alhsin.F, 

asetimj)  NOTES:  (if  cbar<0,  value  computed  in 
coamm.F )  (because  of  numerical  techniques, 
cbar  >  2*IV1) 


ID  array:  Climatological  values  of  relative  humidity  at 
user  specified  pressure  (pr)  levels  (getng.F). 


Speed  of  gravity  waves  cstar  =  {zr/7!:)-^j{g/9^)dd^  jdz 

where  Zt  the  is  domain  top;  used  for  radiative  lateral 
boundaries  (ibdya=jbdya=  4  or  5)  (bwind.F) 


Minimum  horizontal  resolution  for  Kain-Fritsch 
convective  parameterization  scheme  (amodeLF) 


Large,  coarse  mesh  time  step;  used  for  integration  of 
slow  modes  NOTE:  (cfl  condition:  delta  <  dsigma/w 
&  delta  <  0.7*delx/(2*IVI))  (inner  nest  time  steps 


VAR  DEFAULT  VALUE _ DEFINITION _ 

_ computed  internally) 


delx 

40000.0  m  (specify  in 
degrees  of  longitude  when 
nproi=5) 

Coarse  mesh  grid  spacing  in  the  x-direction;  used  for 
computing  horizontal  spatial  derivatives  NOTE: 

(inner  nest  spacing  computed  internally) 

dely 

40000.0  m  (specify  in 
degrees  of  latitude  when 
nproj=5) 

Coarse  mesh  grid  spacing  in  the  y-direction;  used  for 
computing  horizontal  spatial  derivatives  NOTE: 

(inner  nest  spacing  computed  internally) 

diOnd 

0.01 

2nd  order  dimensionless  diffusion  coefficient 
dif2nd=Kd2At/Ax^  <  1/8  for  stability;  used  to  control 
small  scale  numerical  noise  {asetup.F) 

dif4th 

-0.0025 

4th  order  dimensionless  diffusion  coefficient 
dif4th=Kd4At/Ax'^  >  -1/32  for  stability;  used  to  control 
small  scale  numerical  noise  {asetup.F)  NOTE:  (if 
non-  negative,  2nd  order  diffusion  is  used) 

divdmp 

10.0 

Dimensionless  divergence  damping  coefficient 
divdmp=KdivA'i/Ax^  <  1/2  for  stability;  used  to  control 
amplification  of  the  acoustic  mode  (alhs.F)  NOTE: 
(smaller  values  produce  greater  damping) 

dsigma 

500.0  m 

ID  array:  Vertical  grid  spacing  in  the  atmosphere;  used 
to  compute  vertical  spatial  derivatives  NOTE;  (non 
uniform  vertical  spacing  is  permitted) 

dtrad 

3600.0  s 

Radiation  time  interval;  used  to  determine  frequency  of 
radiation  computations  {coamm.F) 

dxmeso 

9000.0  m 

Upper  limit  of  horizontal  resolution  below  which  the 
vertical  advection  of  tke,  raindrops  and  snow  are 
included  in  prognostic  equations  {qfore.F,  qforqx.F) 

dzqlim 

500.0,250.0,100.0,100.0, 
100.0, 100.,  100  m 

Upper  limit  of  height,  for  each  mesh,  below  which 
vertical  spacing  for  vertical  advection  of  rain  and  snow 
is  computed  {dzqrs.F) 

flat 

0.0° 

Latitude;  Used  for  Cartesian  grids  (nproj=4)  to 
compute  the  Coriolis  parameter  {grid.F) 

1.0 

Maximum  ground  wetness  (sfcpar.f) 

hgtsnd 

999.0  m 

ID  array:  Heights  above  surface  of  idealized  sounding 
data  {refsnd.F)  NOTE:  (computed  internally) 

hgtvis 

10.,  50.,  100.,  200. 

,  300.,  400.,  500.,  750. 

,  1000.,1250.,  1500.,  2000. 

,  2500.,3000.,  3500.,  4000. 

,  4500.,5000.,  6000.,  7000. 

,  8000.,9000.,10000.,12000. 
,15000., 75*0.0  m 

ID  array:  Heights  above  the  surface  that  3D  graphics 
data  are  interpolated  to;  used  for  Visual  output  when 
Iprvis=.false.  {visout.F) 

phntl 

30.0° 

Latitude  where  map  projection  intersects  the  earth  and 
map  factor  is  1.0;  used  to  setup  map  projection  (grid.F) 
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VAR 


phnt2 


pr 


prsvis 


psfcO 


psnd 


qsnd 


rdtime 


rlatl 


riongl 


DEFAULT  VALUE 

DEFINITION 

NOTE:  (varies  from  -90°at  the  south  pole  to  +90°  at 
the  north  pole) 

60.0° 

2nd  Latitude  where  map  projection  intersects  the  earth; 
used  to  setup  lambert  conformal  map 
projection(grid.F) 

NOTE:  (varies  from  -90°at  the  south  pole  to  +90°  at 
the  north  pole) 

10.0,  20.0,  30.0,  50.0 
,  70.0,  100.0,150.0,  200.0 
,250.0,  300.0,400.0,  500.0 
,700.0,  850.0,925.0,1000.0 
,  84*0.0  mb 

ID  array:  Pressure  levels  corresponding  to  analyzed 
fields:  used  to  obtain  and  generate  real  data  initial 
conditions 

{getng.F,  getfgl.F,  inlvlF) 

0.0 

Computed  internally 

1030.,1000.,975.,950.,925. 

,  900.,  875.,850.,825.,800. 

,  775.,  750.,725.,700.,650. 

,  600.,  550.,500.,400.,350. 
,300.,  250.,200.,100. 

,76*0.0  mb 

ID  array:  Pressure  levels  that  3D  graphics  data  are 
interpolated  to;  used  for  Visual  output  when 
lprvis=.true.  (visout.F) 

100000.0  kg/ms 

Surface  pressure;  used  to  compute  terrain  pressure  in 
idealized  cases  (user_sfc.F) 

-999.0 

ID  array:  Pressure  or  height  levels  of  idealized 
sounding  data  (refsnd.F)  NOTE:  (pressure  in  mb 
when  ipsnd=l,  height  in  meters  when  ipsnd=2) 

-999.0 

ID  array:  Moisture  of  idealized  sounding  data 
(refsnd.F)  NOTE:  (moisture  in  form  of  dewpoint 
temperature  degrees  Celcius  when  iqsnd=l,  dewpoint 
temperature  degrees  Kelvin  when  iqsnd=2,  mixing 
ratio  (g/kg)  when  iqsnd=3,  relative  humidity  (%)  when 
iqsnd=4) 

120.0  s 

Time  scale  for  damping  u,v,w,0  in  sponge  layer;  used 
to  control  vertical  reflections  {coamm.F)  NOTES: 

(the  larger  the  time  scale,  the  weaker  the  damping) 
(rdtime  =  l/(35  w(2;r/L^ ))  where  Lz  is  the  vertical 
wavelength) 

40.0° 

Reference  latitude,  corresponds  to  reference  point 
(iref(l)Jref(l));  used  to  setup  the  coarse  grid  {grid.F) 
NOTE:  (varies  from  -90°at  the  south  pole  to  +90°  at 
the  north  pole) 

240.0° 

Reference  longitude,  corresponds  to  reference  point 
(iref(l)jjref(l));  used  to  setup  the  coarse  grid  {grid.F) 
NOTE:  (varies  from  -90°at  the  south  pole  to  +90°  at 
the  north  pole) 
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VAR 

DEFAULT  VALUE  I 

it':-'  i.' '  V: 

robert 

0.05 

Coefficient  for  coupling  the  three  time  levels;  used  to 
filter  high  frequency  oscillations  (amodeLF,  atimef.F) 
NOTE:  (suggested  value  of  robert=0.05  when 
Iddamp=.true.,  otherwise,  robert=0. 1-0.2) 

seatmp 

285.0  K 

Idealized  surface  temperature;  used  to  specify  the 
surface  temperature  when  iseaflg  =  0  (user_sfc.F)  ' 

sdmiilt 

0.0 

Multiplication  factor  for  standard  deviation  of  terrain 
height  {sfcpar.F) 

sfcwet 

0.1 

Idealized  value  for  ground  wetness;  used  to  specify  the 
ground  wetness  when  igwtflg  =  o  (user_sfc.F)  NOTE; 
(0.0  =  saturated,  1.0  =  dry) 

shO 

0.675 

TKE  constant;  used  to  compute  the  TKE  coefficient  sh 
when  iashsm=2,3,4  (afore.F)  NOTE:  (typically 
sh0=0.675  when  Ax»Az;  sh0=0.3  when  Ax-Az) 

smO 

0.5 

TKE  constant;  used  to  compute  the  TICE  coefficient 
sm  when  iashsm=2,3,4  (afore.F)  NOTE:  (typically 
sm0=0.5  when  Ax»Az;  sm0=0.1  when  AxsAz) 

1 

Factor  for  stretching  mesh  horizontally  (grid.F) 

tamp 

1.5 

Amplitude  of  idealized  temperature  perturbation 
(fldsOS.F) 

tld 

3600.0  s 

ID  graphics  time  interval;  used  to  determine  frequency 
of  coamld  output  (asavld.F) 

tvisl 

0.0  s 

ID  array:  3D  graphics  time  interval  for  each  mesh; 
used  to  determine  frequency  of  Visual  set  1  output 
(visout.F) 

tvis2 

0.0  s 

ID  array:  3D  graphics  time  interval  for  each  mesh; 
used  to  determine  frequency  of  Visual  set  2  output 
(visout.F) 

tsnd 

-999.0 

ID  array:  Temperatures  of  idealized  sounding  data 
(refsnd.F)  NOTE:  (temperature  in  degrees  Kelvin 
when  itsnd=l,  in  degrees  Celcius  when  itsnd=2,  in 
theta  when  itsnd=3) 

ugeoa 

vgeoa 

0.0  m/s 

ID  arrays:  Geostrophic  wind  components;  used  to 
compute  the  synoptic  scale  pressure  gradient  force 
when  lvgeo=.true.  (arhsu.F,  arhsv.F) 

umean 

vmean 

0.0  m/s 

Constant  values  representing  a  mean  wind;  used  to 
subtract  a  mean  wind  from  the  idealized  sounding  data 
(istate.F) 

rmm 

-999.0  m/s 

ID  arrays:  Winds  of  idealized  sounding  data  (refsnd.F) 
NOTE:  (winds  in  terms  of  components  (m/s)  when 
iuvsnd=2,  in  terms  of  direction-usnd  and  speed-vsnd 
when  iuvsnd=l) 

zOzOmx 

0.1  m 

Maximum  surface  roughness  (sfcpar.F) 
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VAR  DEFAULT  VALUE 


zrough  0.01  m 


INTEGER  VARIABLES 


iadvct  1 


ialbflg  1 


iahsgm  1 


iamxgl  1 


iashsm  1 


ibdya  1 


lease  I  1 


DEFINITION 


Idealized  value  for  surface  roughness;  used  to  specify 
the  surface  roughness  when  izOflg  =  0  (user_sfc.F) 
NOTE:  (0.0  =  smooth,  1.0  =  rough) 


Advection  type  identifier  {arhs*.F,  afor*.F,  advsc.F): 

(1)  4th  order  advective 

(2)  2nd  order  advective 

(3)  2nd  order  advective  with  2nd  order  upstream  for  0 

(4)  spectral 

(5)  2nd  order  flux 

(6)  mixed:  2nd  order  advective  and  flux 


Indicator  for  use  of  idealized  albedo  data  for  each  mesh 
(user_rfc.F)  (NOTE:  albedo  given  by  idealized  value 
(albdo)  when  ialbflg=0,  otherwise  uses  database) 


Horizontal  subgrid  scale  mixing  identifier: 

(arhs*.F,  afor*.F,  hmix*.F) 

(1)  eddy  mixing  coefficient  based  upon  deformation 
field  (use  when  Ax»Az) 

(2)  eddy  mixing  coefficient  based  upon  TKE 
prediction  (use  when  AxsAz) 


Mixing  length  identifier:  options  for  the  mixing 
lengths  all,  al2,  aid  are  given  in  Table  D-1  below 
(afore.F,coamm.F) 


TKE  coefficients  sm  and  sh  identifier;  used  to  define 
eddy  mixing  coefficients  (afore.F): 

(1)  Ri  dependent  (Mellor  &  Yamada,  74) 

(2)  sin=sm0,  sh=sh0 

(3)  sm=sm0,  sh=(  1 .0+2.0  al2/all)*sm0 

(4) sm=sm0,  sh=2.13*sm0 

(5)  Ri  dependent  (Mellor  &  Yamada,  82) 


x-direction  boundary  condition  identifier  (abdy*.F): 
idealized 

(1)  fixed 

(2)  extrapolated 

(3)  periodic 

(4)  radiation 

(5)  radiation  and  extrapolation 

real  data 

(6)  Perkey-Krietzberg 

(7)  Davies 


Case  number  identifier;  used  to  indicate  real  or 
idealized  da.ta  (coama.F,  flds##.F): 

(0)  real  data 

(##)  value  corresponds  to  the  number  in  flds##.F 


VAR 

DEFAULT  VALUE 

DEFCTOON  ■ 

filename,  a  file  where  perturbations  of  the  mean  state 
are  specified 

icycle 

Ohr. 

Frequency  of  a  data  assimilation  cycle;  used  to 
read/write  input/output  fields  (coama.F,  coamm.F) 

idgadj 

2 

x-grid  point;  used  for  printout  of  microphysical 
diagnostics  {adjtq.F) 

ifsave 

1 

Indicator  for  saving  sigma  level  data  (coamm.F): 

(1)  at  specific  times  given  by  ksaves 

(2)  at  time  intervals  given  by  isavfrq 

igwtflg 

1 

Indicator  for  use  of  idealized  ground  wetness  for  each 
mesh  {user_sfc.F)  (NOTE:  ground  wemess  given  by 
idealized  value  (sfcwet)  when  igwtflg=0,  otherwise 
uses  database) 

ii 

21 

ID  array:  Parent  meshes  x-grid  point  corresponding  to 
left  boundary  of  next  inner  nest;  used  to  setup  inner 
nests  (coama.F) 

ilndflg 

1 

Indicator  for  use  of  idealized  land/sea  grid  point  for 
each  mesh  (userjsfc.F)  (NOTE:  land/sea  point  given 
by  idealized  value  (alndpct)  when  ilndflg=0, 
otherwise  uses  database) 

iloi 

0 

MVOI  volume  type  identifier  (coama.F)‘. 

(0)  ivol,  jvol  =  number  of  MVOI  volumes 

(1)  ivol,  jvol  =  distance  in  kilometers  of  an  MVOI 

volume 

ilsp 

1 

Indicator  for  use  of  explicit  moist  physics  for  each 
mesh  (amodeLF)  NOTE:  (explicit  moist  physics 
when  ilsp=l,  large  scale  precip  param  when  ilsp=2) 

ipsnd 

1 

Indicator  of  units  of  idealized  pressure  data  (refsnd.F) 
NOTE:  (pressure  in  mb  when  ipsnd=l,  height  in 
meters  when  ipsnd=2) 

iqsnd 

1 

Indicator  of  units  of  idealized  moisture  data  (refsnd.F) 
NOTE:  (moisture  in  form  of  dewpoint  temperature 
degrees  Celcius  when  iqsnd=l,  dewpoint  temperature 
degrees  Kelvin  when  iqsnd=2,  mixing  ratio  (g/kg) 
when  iqsnd=3,  relative  humidity  (%)  when  iqsnd=4) 

iref 

55 

ID  array:  x-grid  point  on  each  mesh  at  reference 
latitude  rlatl;  used  to  setup  the  model  domains 
(grid.F)  NOTE:  (iref(l)  can  be  larger  than  ma(l) 
when  alnnt  and  rlongl  fall  outside  the  model  domain  - 
-  for  rotated  grids)  (iref  for  inner  nests  computed 
internally) 

isavefrq 

-1,0,0  hr,min,sec 

Times  for  saving  sigma  level  data;  used  to  graphically 
view  CO  AMPS  forecast  fields  (coamm.F)  NOTE: 
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VAR 

DEFAULT  VALUE 

DEFINITION 

(used  when  ifsave=2) 

iseaflg 

1 

Indicator  for  use  of  idealized  surface  temperature  data 
for  each  mesh  {user_sfc.F)  (NOTE:  surface 
temperature  given  by  idealized  value  (seatmp)  when 
iseaflg=0,  otherwise  uses  database) 

istrl 

1 

Factor  for  horizontally  stretched  grid  {grid.F) 

istr2 

1 

Factor  for  horizontally  stretched  grid  (^rid.F) 

itopoflg 

1 

Indicator  for  use  of  idealized  terrain  data  for  each  mesh 
(user_sfc.F)  (NOTE:  topography  height  set  to  zero 
when  itopoflg=0,  otherwise  uses  database) 

itauf 

12  hr. 

Ending  time  for  obtaining  and  writing  out  coarse  mesh 
boundary  tendencies;  used  for  data  assimilation 
(coama.F) 

itauin 

12  hr. 

Initial  time  (and  frequency)  for  obtaining  and  writing 
out  coarse  mesh  boundary  tendencies  from  NOGAPS 
fields;  used  for  data  assimilation  (coama.F)  NOTE: 
(value  should  be  equivalent  to  forecast  interval  of 
NOGAPS  fields) 

itsnd 

1 

Indicator  of  units  of  idealized  temperature  data 
(refsnd.F)  NOTE:  (temperature  in  degrees  Kelvin 
when  itsnd=l,  in  degrees  Celcius  when  itsnd=2,  in 
theta  when  itsnd=3) 

iupd 

0 

Data  assimilation  identifier  {coama.F,  asetup.F, 
inlvl.F,  coamm.F): 

(0)  no  update;  analyzed  NOGAPS  fields  on  p  levels 

(1)  full  update;  analyzed  CO  AMPS  fields  on  p  levels 

(2)  incremental  update;  CO  AMPS  fields  on  sigma-z 
levels  with  MVOI  increments  interpolated  to  sigma-z 
levels 

iuvsnd 

1 

Indicator  of  units  of  idealized  wind  data  {refsnd.F) 
NOTE;  (winds  in  terms  of  components  (m/s)  when 
iuvsnd=2,  in  terms  of  direction-usnd  and  speed-vsnd 
when  iuvsnd=l) 

ivar 

0 

Nonlinear  variational  balancing  identifier;  used  in 

MVOI  analysis  (oi_anl.F): 

(0)  no  balancing 
(1)  balancing 

ivol 

5 

ID  array;  Number  of  regions  (MVOI  volumes)  that 
each  mesh  is  divided  into  in  the  x-direction;  used  for 
MVOI  analysis  {oivol.F,  oi_anl.F) 

ild 

0 

2D  array:  ID  graphics  x-grid  points  for  each  mesh; 
used  to  specify  location  of  ID  graphics  output 
{asavld.F) 
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VAR  1  DEEAXILT  VALUE  ..  I  IHMmm 


izOflg 

1 

Indicator  for  use  of  idealized  surface  roughness  data 
for  each  mesh  {userjsfc.F)  (NOTE:  surface  roughness 
set  to  idealized  value  (zrough)  when  iz0flg=0, 
otherwise  uses  database) 

jbdya 

1 

y-direction  boundary  condition  identifier  (abdy*F): 
idealized 

(1)  fixed 

(2)  extrapolated 

(3)  periodic 

(4)  radiation 

(5)  radiation  and  extrapolation 

real  data 

(6)  Perkey-Krietzberg 

(7)  Davies 

NOTE:  (in  general  ibdya=jbdya) 

jdgadj 

2 

y-grid  point:  used  for  printout  of  microphysical 
diagnostics  (adjtq.F) 

jref 

42 

ID  array:  y-grid  point  at  reference  longitude  riongl; 
used  to  setup  the  model  domains  (grid.F) 

NOTE:  (jref  for  inner  nests  computed  internally) 

jj 

21 

ID  array:  Parent  meshes  y-grid  point  corresponding  to 
lower  boundary  of  next  inner  nest;  used  to  setup  inner 
nests  {coama.F) 

jqq 

550 

Max  number  of  observations  per  MVOI  volume 
{oi_anl.F) 

jvol 

5 

ID  array:  Number  of  regions  (MVOI  volumes)  that 
each  mesh  is  divided  into  in  the  y-direction;  used  for 
MVOI  analysis  {oivolF,  oi_anl.F) 

jld 

0 

2D  array:  ID  graphics  y-grid  points  for  each  mesh; 
used  to  specify  location  of  ID  graphics  output 
(asavld.F) 

kgetbc 

12  hr. 

Frequency  for  updating  coarse  mesh  boundary 
tendencies;  used  in  data  assimilation  {coamm.F) 

NOTE:  (value  chosen  so  that  kgetbc/itauin  is  a 
whole  number) 

kprnta 

6,0,0  hr,min,sec 

Frequency  for  writing  forecast  fields  to  standard 
output;  used  for  viewing  slices  of  numerical  model 
data  {coamm.F) 

ksavea 

-1,0,0  hr,min,sec 

Frequency  for  writing  restart  files  to  directory  dsnrff; 
used  to  obtain  fields  used  for  model  restart  at  a 
nonzero  time  {coamm.F) 

ksaves 

-1  lir,min,sec 

2D  array:  Times  for  saving  sigma  level  data  for  each 
mesh;  used  to  graphically  view  CO  AMPS  forecast 
fields  {coamm.F)  NOTE:  (used  when  ifsave=l) 
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VAR 

DEFAULT  VALUE 

DEFINITION 

ktauf 

12,0,0  hr,min,sec 

2D  array:  Ending  time  of  forecast  for  each  mesh 
{coamm.F) 

ktaust 

0,0,0  hr,min,sec 

Beginning  time  of  forecast  (coainm.f)  NOTE:  (ktaust 
usually  equal  to  zero  except  for  a  model  restart  -  see 
ksavea  to  save  fields  for  a  model  restart) 

mtaua 

5 

Number  of  small  coarse  mesh  time  steps  per  (delta) 
large  time  steps;  used  for  integration  of  fast  modes 
(alhs.F)  NOTE:  (cfl  condition:  deltaua  <  dsigma/w 
&  deltaua  <  0.7*delx/(cbar+IVI) )  (inner  nest  time 
steps  computed  internally)  (because  of  numerical 
techniques  used,  mtaua  typically  set  to  2) 

nftopo 

3 

Number  of  passes  through  25-point  filter;  used  to 
smooth  terrain  field  (sfcpar.F) 

nfpass 

3 

Number  of  passes  through  9-  or  25-point  filter;  used  to 
smooth  output  and  analysis  fields  (filt9.F ,  filtlS.F) 

njump 

3 

njump- 1  is  number  of  x-  and  y-grid  points  skipped  in 
radiation  computations  (radiat.F) 

noisub 

12 

Number  of  processors;  for  MVOI  multitasking 
(oi_anl.F) 

nproj 

2 

Map  projection  identifier;  used  to  setup  the  horizontal 
domain  (grid.F): 

(1)  mercator 

(2)  lambert  conformal 

(3)  polar  stereographic 

(4)  Cartesian  coordinates 

(5)  spherical 

nrdamp  8  Number  of  upper  model  levels  that  are  in  the  sponge 

layer;  used  to  damp  vertical  reflections  (coamm.F) 
_ NOTE:  (8  or  more  layers  recommended) _ 


LOGICAL  VARIABLES 


laper 

•false. 

Check  for  aperiodicity  at  boundaries;  used  as  a  check 
for  idealized  cases  with  periodic  boundary  conditions 
(ibdya=jbdya=3)  {asyml.F,  asym2.F): 
true  -  check  periodicity 
false  -  do  not  do  above 

lasym 

•false. 

Check  for  asymmetry  in  fields;  used  as  a  check  for 

idealized  cases  with  symmetric  initial  conditions  and 

forcing  {amodeLF): 

true  -  check  for  asymmetry 

false  -  do  not  do  above 

Icool 

•false. 

Simple  radiational  cooling  {aforth.F): 
true  -  cool  2°  per  day  below  10  km 
false  -  do  not  do  above 
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DEFAULT  VALUE 

DEFINITION 

Icupar 

•false. 

Kain-Fritseh  convective  parameterization  (amodelF): 
true  -  perform  convective  adjustments  to  (0,  qv,  qc,  qr) 
on  meshes  with  resolution  Ax>cupdx 
false  -  do  not  do  above 

Icuppr 

•false. 

Kain-Fritsch  parameterization  diagnostics  {kfpam.F): 
true  -  print  diagnostics 
false  -  do  not  do  above 

Idbms 

•false. 

Read  NOGAPS  initial  conditions  (coama.F)\ 
true  -  ISIS  database 
false  -  flat  files 

Idbmsi 

•false. 

Read/write  CO  AMPS  input/output  {coama.F, 

coamm.F): 

true  -  ISIS  database 

false  -  flat  files 

Idbmso 

•false. 

Write  CO  AMPS  analysis  fields  {coama.F)\ 
true  -  ISIS  database 
false  -  flat  files 

Iddamp 

•false. 

Divergence  damping;  used  to  reduce  accumulation  of 
mass  within  a  grid  volume  (alhs.F): 
true  -  use  divergence  damping 
false  -  do  not  do  above 

Idgadj 

•false. 

Microphysical  diagnostics  (adjtq.F): 
true  -  print  diagnostics 
false  -  do  not  do  above 

Idiff 

•true. 

Numerical  diffusion  (arhs*.F,  afor*F): 
true  -  use  diffusion 
false  -  do  not  do  above 

Ifftf 

•false. 

FFT  filter  —  spectral  truncation  of  wavelengths  <  3Ax; 
used  in  place  of  numerical  diffusion  for  LES  cases 
with  periodic  boundary  conditions  (ibdya=jbdya=3) 

mn-F)-. 

true  -  use  FFT  filter  (requires  that  Idiff  =.false.  also) 
false  -  do  not  do  above 

Iflux 

•true. 

Surface  fluxes  {sfcflx.F): 

true  -  nonzero  surface  heat,  moisture,  momentum  flux 
false  -  zero  surface  heat,  moisture,  momentum  flux 

Ues 

•false. 

Large  eddy  simulation;  use  when  Ax=Az  {afore.F, 
sfcflx.F): 

true  -  use  terms  and  coefficients  in  the  TKE  and 

surface  flux  computations  that  are  appropriate 
for  3D  turbulence 
false  -  do  not  do  above 

lice 

•true. 

Ice  physics  {adjtq.F,  amodeLF): 

D-11 


VAR  DEFAULT  VALUE _ DEFTNITION _ 

true  -  temporal  variations  and  adjustments  computed 
for  frozen  moisture  variables  (qi,  qs)  (requires 
that  lmoist=.true.  also) 

_ false  -  set  frozen  moisture  variables  to  zero _ 

Unit  .false.  Initialize  certain  variables  {ainit.F): 

Real  data 

true  -  eliminate  net  inflow/outflow,  compute  vertical 
velocity 

false  -  do  not  do  above 
Idealized  data  (based  upon  value  of  icase) 
true  -  initialize  pressure,  or  potential  temperature,  or 
add  perturbations  to  horizontal  winds 

_ false  -  do  not  do  above _ 

Imoist  .true.  Moist  physics  (adjtq.F,  amodelF): 

true  -  temporal  variations  and  adjustments  computed 
for  liquid  moisture  variables  (qc,  qr)  and  frozen 
moisture  variables  (qj,  qs)  when  Uce=.true.  also) 

_ false  -  set  liquid  and  frozen  moisture  variables  to  zero 

Inrhty  .true.  Initialization  of  surface  fields  (gethty.F,  nethty.F): 

true  -  use  CO  AMPS  history  fields  if  available 

_ false  -  use  other  source:  NOGAPS,  climatology  fields 

locard  .false.  Read  in  2D  horizontal  level  plotting  information 

{coama.F,  coamm.F): 

true  -  ISIS  operational  format 

_ false  -  do  not  do  above _ 

loi  .false.  Multivariate  optimum  interpolation  (MVOI)  analysis 

(oi_anl.F,  anlfld.F): 
true  -  perform  MVOI 

_ false  -  do  not  do  above _ 

loimf  .false.  Inner  nest  MVOI  analysis  (o/_an/.F,  an/^<i.F): 

true  -  perform  MVOI  on  inner  nest  fields 
false  -  interpolate  coarse  mesh  increments  to  inner 

_ nest  fields _ 

Iprint  .false.  Surface  parameter  diagnostics  (sfcpar.F): 

true  -  print  diagnostics 

_ false  -  do  not  do  above _ 

Iprvis  .false.  Visual  data  vertical  interpolation  (yisout.F): 

true  -  pressure  levels  given  by  the  prsvis  array 
false  -  height  levels  (above  surface)  given  by  the 

_ hgtvis  array _ 

Iqanl  .false.  Moisture  analysis  {anlfld.F,  oi_anl.F): 

true  -  perform  analysis  of  initial  moisture  field 
false  -  do  not  do  above 


VAR 

DEFAULT  VAIAIE;  , 

Irad 

.false. 

Long  and  short  wave  atmospheric  radiation  scheme 
based  on  Harshvardhan  {coamm.F,  amodelF): 
true  -  radiative  adjustments  to  0  computed  when 

integration  loop  coincides  with  a  dtrad  time 
interval 

false  -  do  not  do  above 

Iralee 

.false. 

Raleigh  damping  layer  {coamm.F): 
true  -  damp  (u,v,w,0)  to  mean  state  in  the  top 
nrdamp  levels 
false  -  do  not  do  above 

Isfcen 

.false. 

Surface  energy  budget  (sfcenr.F): 

true  -  temporal  changes  in  ground  temperature  and 

wetness  based  upon  surface  fluxes  and  radiation 
false  -  no  temporal  changes  in  ground  temperature  or 
wetness 

Ispong 

.false. 

Sponge  layer  (coamm.F): 

true  -  damp  (u,v,w,0)  to  smoothed  values  in  the  top 
nrdamp  levels 
false  -  do  not  do  above 

Itanl 

.false. 

Temperature  analysis  {anlfld.F,  oi_anl.F): 
true  -  perform  analysis  of  initial  temperature  field 
false  -  do  not  do  above 

Itanls 

.false. 

Surface  temperature  analysis  {anlfld.F,  oi_anl.F): 
true  -  perform  analysis  of  initial  surface  temperature 
field 

false  -  do  not  do  above 

Itke 

.true. 

Subgrid  scale  mixing  {afore.F): 

true  -  temporal  variations  in  TKE  computed 

false  -  no  temporal  variations  in  TKE 

Itopoa 

.true. 

Topography  {coama.F,  insfc.F): 

true  -  terrain  field  given  by  user  or  by  data  base 

false  -terrain  field  set  to  zero 

Ivarad 

.false. 

Hydrostatic  variational  adjustment  {stdp2z.F, 
incrup.F): 

true  -  adjust  pressure  and  potential  temperature 
perturbations 
false  -  do  not  do  above 

Ivgeo 

.false. 

Geostrophic  winds  {arhsu.F,  arhsv.F): 
true  -  include  large  scale  pressure  gradient  force 
computed  from  ugeo  and  vgeo  arrays 
false  -  do  not  do  above 

VAR 

DEFAULT  VALUE 

DEFINITION 

Ivpert 

•false. 

Perturbation  added  to  v-component  (ainit.F): 
true  -  add  perturbations 
false  -  do  not  do  above 

Iwritu 

•false. 

Formatting  of  model  input/output  data: 
true  -  unformatted  (32-bit  IEEE) 
false  -  formatted  (5el3.6) 

Table  D-1 :  Namelist  input  options  for  the  TKE  mixing  length  parameter  iamxgl.  Column  one  contains  the 


variables  and  arrays  in  the  TKE  subroutine  afore.F  that  are  controlled  by  specifying  values  (l)-(5) 
for  iamxgl.  The  alO  is  a  user  defined  input  and  alOO  depends  on  the  zero  and  first  mornents  of  TKE. 


al2= 

(diagnostically 
computed 
mixing  length) 

7fj  same  aI2  as  (1)  same  aI2  as  (1)  fna(blht,  stab,  TKE)  alOO  •  zK 

755751  (a;00+z«r) 

all  -  (unstable) 

all= 

(used  in  al2; 
constant  wrt 
time) 

alQ  •  ZK  (AxAyAz)^^^  alO*Az  all  not  used  all  not  used 

{alQ+zK) 

ald= 

(used  in  the  TKE 
dissipation  term) 

al2  same  aid  as  (1)  same  aid  as  (1)  fnb(blht,  stab,  TKE)  al2 

coef= 

(used  in  the  TKE 
dissipation  term) 

0  19+0  74 

all  (1) 

reference 

Mellor  &  Yamada  Klemp  &  —  Therry  &  Lacarrere  Mellor  & 

1974  Wilhelmson  1983  Yamada  1982 

1978 

App6ndix  E:  Three  sample  map  projections,  input  parameters,  and  model  domains 


10*  20*  30*  40 

LATITUDE 


70*  80*  90* 


Figure  E-1 :  Schematic  representation  of  (a)  a  mercator  map  projection  with  standard 
latitude  at  20°N  (Haltiner  and  Williams,  1980),  and  (b)  the  corresponding  graph 
showing  variations  in  map  scale  factor  with  latitude  (Saucier,  1955). 


Table  E-1:  Namelist  input  parameters  and  values  for  the  mercator 
_ map  projection  shown  in  Figure  E-2. 


Namelist  Parameters  I  Definition  and  Use 


not  used 


coarse  mesh  (1)  arid  soacina  in  the  x-direction 


coarse  mesh  (1)  arid  spacinq  in  the  v-direction 


coarse  mesh  (1)  x-grid  point  corresponding  to  the  reference  latitude 


coarse  mesh  (1)  x-grid  point  corresponding  to  left  boundary  of  inner  nest  (2 


not  used 


rid  point  corresponding  to  the  reference  longitude 


rid  point  corresponding  to  left  boundary  of  inner  nest  (2 


not  used 


coarse  mesh  (1)  x-grid  points 


inner  nest  (2)  x-grid  points 


inner  nest  (3)  x-grid  points 


rid  points 


rid  points 


rid  points 


model  domains 


mecator  map  projection 


deix  =  120000.0, 


IWMS 


nnest  = 


3 

0. 

latitude  where  ma 

P  P 

rejection  intersects  the  earth  (i.e.  map  factor=1 

not  used 


reference  latitude 


reference  longitude 


E-2 


Figure  E-2:  Example  of  a  COAMPS  horizontal  domain  using  a  mercator  map  projection  over 
the  western  Pacific,  The  arrows  correlate  locations  on  the  grid  to  the  input  parameters  and  their 
user-specified  value.  These  pararheters  are  further  defined  in  Table  E-1, 


Figure  E-3:  Schematic  representation  of  (a)  a  lambert  conformal  map  projection  with 
standard  latitudes  at  30°N  and  60°N  (Haitiner  and  Williams,  1980),  and  (b)  the 
corresponding  graph  showing  variations  in  map  scale  factor  with  latitude  (Saucier,  1955). 


t 


Table  E 

-2;  Namelist  input  parameters  and  values  for  the  lambert 
conformal  map  proi action  shown  in  Figure  E-4. 

Namelist  Parameters 

Definition  and  Use 

alnnt 

tz 

16.5, 

north/south  longitude  used  to  orient  and  align  the  grid 

deix 

= 

81000.0, 

coarse  mesh  (1)  grid  spacing  in  the  x-direction 

dely 

z= 

81000.0, 

coarse  mesh  (1)  grid  spacing  in  the  y-direction 

irefd) 

=z 

31, 

coarse  mesh  (1)  x-qrid  point  corresponding  to  the  reference  latitude 

ii(2) 

= 

22, 

coarse  mesh  (1)  x-grid  point  corresponding  to  left  boundary  of  inner  nest  (2) 

ii{3) 

=z 

16, 

inner  nest  (2)  x-grid  point  corresponding  to  left  boundary  of  inner  nest  (3) 

irefd) 

= 

31, 

coarse  mesh  (1)  v-grid  point  corresponding  to  the  reference  longitude 

ii(2) 

= 

22, 

coarse  mesh  (1)  y-grid  point  corresponding  to  left  boundary  of  inner  nest  (2) 

ii(3) 

= 

16, 

inner  nest  (2)  y-grid  point  corresponding  to  left  boundary  of  inner  nest  (3) 

ma(1) 

= 

61, 

coarse  mesh  (1)  x-qrid  points 

ma(2) 

= 

61, 

inner  nest  (2)  x-grid  points 

ma(3) 

= 

91, 

inner  nest  (3)  x-grid  points 

nad) 

61, 

coarse  mesh  (1)  y-grid  points 

na(2) 

= 

61, 

inner  nest  (2)  v-grid  points 

na(3) 

= 

91, 

inner  nest  (3)  y-grid  points 

nnest 

3, 

model  domains 

nproj 

= 

2, 

lambert  conformal  map  projection 

phnt1 

= 

60.0, 

latitude  where  map  projection  intersects  the  earth  (i.e.  map  factor=1) 

phnt2 

30.0, 

second  latitude  where  map  projection  intersects  the  earth 

riatl 

= 

42.5, 

reference  latitude 

im 

= 

16.5, 

reference  longitude 
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Figure  E-5:  Schematic  representation  of  (a)  a  polar  stereographic  map  projection  with 
standard  latitude  at  90°N  (Haltiner  and  Williams,  1980),  and  (b)  the  corresponding 
graph  showing  variations  in  map  scale  factor  with  latitude  (Saucier,  1955). 


Table  E-3 :  Namelist  input  parameters  and  values  for  the  polar 
_ stereographic  map  projection  shown  in  Figure  E-6. 


Namelist  Parameters 

Definition  and  Use 

ainnt 

280.0, 

north/south  longitude  used  to  orient  and  align  the  grid 

deix 

381000.0, 

coarse  mesh  (1)  grid  spacing  in  the  x-direction 

dely 

= 

coarse  mesh  (1)  grid  spacing  in  the  y-direction 

ihl 

32, 

coarse  mesh  (1)  x-grid  point  corresponding  to  the  standard  latitude 

ii21 

11. 

coarse  mesh  (1)  x-grid  point  corresponding  to  left  boundary  of  inner  nest  (2) 

ii32 

= 

21, 

inner  nest  (2)  x-grid  point  corresponding  to  left  boundary  of  inner  nest  (3) 

ih1 

= 

32, 

coarse  mesh  (1)  y-grid  point  corresponding  to  the  standard  longitude 

ii21 

zz 

11. 

coarse  mesh  (1)  y-grid  point  corresponding  to  left  boundary  of  inner  nest  (2) 

ii32 

= 

21. 

inner  nest  (2)  y-grid  point  corresponding  to  left  boundary  of  inner  nest  (3) 

mal 

=r 

61, 

coarse  mesh  (1)  x-grid  points 

ma2 

= 

55, 

inner  nest  (2)  x-grid  points 

ma3 

=r 

55, 

inner  nest  (3)  x-grid  points 

na1 

= 

61, 

na2 

55. 

na3 

=: 

55, 

nnest 

= 

3, 

model  domains  I 

nproj 

zr 

3, 

laassi 

= 

60.0, 

latitude  where  map  projection  intersects  the  earth  (i.e.  map  factor=1) 

zz 

— —  j 

not  used 

riatl 

zz 

90.0, 

standard  latitude 

IKU 

=r 

280.0, 

standard  longitude 

E-5 


E-6 


Figure  E-6:  Example  of  a  COAMPS  horizontal  domain  using  a  polar  stereographic  map  projection 
over  the  North  Pole.  The  arrows  correlate  locations  on  the  grid  to  the  input  parameters  and  their 
user-specified  value.  These  parameters  are  further  defined  in  Table  E-3. 


